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Abstract 

We consider existence and stability properties of nonlinear spatially periodic or 
quasiperiodic standing waves (SWs) in one-dimensional lattices of coupled anhar- 
monic oscillators. Specifically, we consider Klein-Gordon (KG) chains with either 
soft (e.g., Morse) or hard (e.g., quartic) on-site potentials, as well as discrete non- 
linear Schrodinger (DNLS) chains approximating the small-amplitude dynamics of 
KG chains with weak inter-site coupling. The SWs are constructed as exact time- 
periodic multibreather solutions from the anticontinuous limit of uncoupled oscil- 
lators. In the validity regime of the DNLS approximation these solutions can be 
continued into the linear phonon band, where they merge into standard harmonic 
SWs. For SWs with incommensurate wave vectors, this continuation is associated 
with an inverse transition by breaking of analyticity. When the DNLS approximation 
is not valid, the continuation may be interrupted by bifurcations associated with 
resonances with higher harmonics of the SW. Concerning the stability, we identify 
one class of SWs which are always linearly stable close to the anticontinuous limit. 
However, approaching the linear limit all SWs with nontrivial wave vectors become 
unstable through oscillatory instabilities, persisting for arbitrarily small amplitudes 
in infinite lattices. Investigating the dynamics resulting from these instabilities, we 
find two qualitatively different regimes for wave vectors smaller than or larger than 
7r/2, respectively. In one regime persisting breathers are found, while in the other 
regime the system rapidly thermalizes. 
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1 Introduction 



It is well-known that the dynamics of spatially periodic harmonic systems is 
non-ergodic, since it can be completely described as a linear combination (with 
time-independent coefficients) of travelling waves of the form cos{Qx — uj{Q)t) 
with wave vector Q and frequency uj{Q). When there are (even small) anhar- 
monic terms present, most harmonic solutions with amplitude zero are very 
unlikely to be continuable to finite amplitude, because it is generally believed 
(despite the famous work of Fermi, Pasta and Ulam [1]) that anharmonic 
terms in general introduce some ergodicity in the system, which in the long 
term should lead to its thermalization. 

It is nevertheless an interesting question to understand whether special har- 
monic solutions could persist as exact solutions when the system becomes 
anharmonic. Actually, it has been proven recently that small-amplitude peri- 
odic travelling waves (with the form cos{Qx — ujt) in the harmonic limit) do 
exist for Hamiltonian chains of linearly coupled nonlinear oscillators {Klein- 
Gordon (KG) lattices) [2] as well as for chains of mass particles with nonlinear 
nearest-neighbour interactions {Fermi- Pasta- Ulam (FPU) lattices) [3]. How- 
ever, for KG lattices the unique continuation of harmonic travelling waves is 
possible only in a regime of weak inter-site coupling (i.e., close enough to the 
so-called anticontinuous limit [4] of uncoupled oscillators); beyond this regime, 
the situation is more complex and there are many bifurcations where new so- 
lutions appear (such as nanopterons, which are travelling localized waves with 
a noncancelling periodic tail at infinity) [2,5]. In the simpler regime close to 
the anticontinuous limit, numerical techniques for finding time-periodic prop- 
agating waves, based on the concept of discrete breathers (see e.g. [4,6] for 
reviews of recent developments on discrete breathers), were implemented in 
[7]. The solutions were constructed as multi-site breathers {^multibreathers^) 
from the trivial solution at the anticontinuous limit where all oscillators are 
oscillating in phase with the same frequency. An arbitrary spatially uniform 
phase torsion was then applied, and a Newton scheme used to continue these 
solutions to non-zero coupling according to the theory described in [4]. 

However, as was found already in [8] from a discrete nonlinear Schrodinger 
(DNLS) approximation (generally valid for small-amplitude oscillations and 
small couphng, see also Sec. 2.2 below), propagating waves in KG lattices may 
exhibit a modulational (Benjamin-Feir) instability by which an initially uni- 
form wave breaks up into an array of localized pulses. Depending on whether 
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the KG on-site potential is soft or hard (i.e. whether the oscillation frequency 
decreases or increases with increasing oscillation amplitude for the individual 
oscillators), the modulational instability occurs for either small or large wave 
vectors of the travelling wave. (A similar result has been found also for FPU 
lattices [9].) This instability also appears a consequence of the theory of effec- 
tive action for multibreathers with phase torsion as developed in [4] (see also 
further developments in [10]). 

In this paper, we are interested in the continuation of standing waves (SWs), 
which in the harmonic limit are superpositions of two travelling waves with the 
same amplitude and opposite wave vectors, and thus have the form 
cos(Qa;) cos{ujt). We shall here study such solutions for discrete KG systems, 
and for the related DNLS system describing small-amplitude solutions close 
to the anticontinuous limit. Our approach for finding these solutions is simi- 
lar to that used in [11] for disordered KG lattices: we construct appropriate 
multibrcather solutions with frequency outside the linear phonon band at the 
anticontinuous limit and continue them to nonzero coupling; these solutions 
are then continued versus frequency inside the phonon band. We shall pro- 
pose an ansatz for finding the appropriate coding sequence [12,4] for these 
multibreathers, allowing for a smooth continuation of the multibreather ver- 
sus frequency down to zero amplitude as a linear SW. When the spatial period 
of the SW is incommensurate with the lattice period, its envelope (hull) func- 
tion undergoes a 'transition by breaking of analyticity' (TBA) [13-15] at some 
critical value of its frequency. 

Investigating the linear stability of these nonlinear SWs in infinite lattices, we 
find the rather striking result that although there exist linearly stable SWs 
close to the anticontinuous limit, all SWs with nontrivial wave vectors (i.e. 
different from or tt) become linearly unstable through an oscillatory insta- 
bility bifurcation when approaching the linear limit. Thus, small-amplitude 
nonlinear SWs will be unstable also for wave numbers where the correspond- 
ing propagating waves are modulationally stable. These extended instabilities 
are a direct consequence of discreteness, and have not been found within ear- 
lier continuum approximations [16]. However, similarly as found e.g. in [17,18], 
the instabilities are very sensitive to finite size effects, so that in a finite system 
linear stability will generally be recovered close to the linear limit. 

Moreover, we show that when their hull function is non-analytic, these SWs 
can be viewed as arrays of discommensurations (see, e.g., [19]) of a simple 
stationary phonon with wave vector tt or 0, and that each discommensuration 
is unstable (as was found already in [18] for the DNLS model). The effect of 
the instability on a single discommensuration is analysed with the technique 
developed in [20,21], and it is found that, after an initial oscillatory regime, it 
moves spontaneously in the lattice with a rather well-defined velocity, analo- 
gously to the behaviour of the discrete dark mode in the DNLS model found in 
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[22,18]. The long-time behavioTir of the Tinstable SWs can then be interpreted 
as the resuh of inelastic interactions between individual moving discommen- 
surations. As for the DNLS model the SW with wave vector ■n/2 is found 
to exactly coincide with the hne of phase transition recently discovered [23] 
in this model, SWs with wave vectors larger or smaller than 7r/2 belong to 
different phases and therefore yield qualitatively different asymptotic dynam- 
ics: one phase corresponds to normal thermalization, while in the other phase 
persisting large-amplitude localized excitations occur. As we will find from 
numerical simulations, a similar scenario also occurs for KG chains with small 
inter-site coupling. 



The layout of this paper is as follows. In Section 2, we describe the KG model 
and derive a DNLS-type equation for its small- amplitude dynamics. This equa- 
tion generally contains long-range nonlinear interactions, but reduces to the 
ordinary nearest-neighbour DNLS equation in the limit of small coupling. In 
Section 3, we describe the representation of nonlinear SWs as trajectories of 
a symplectic map, which in the DNLS approximation reduces to a cubic map 
in the 2-dimensional real plane [24]. The breakup of Kol'mogorov-Arnol'd- 
Moser (KAM) tori for this map corresponds to a TBA for the hull function 
of incommensurate SWs, and the continuation of periodic and quasiperiodic 
map trajectories to the anticontinuous limit yields the SW multibreather cod- 
ing sequences. In Section 4, we study numerically and analytically the linear 
stability properties of SWs within the DNLS approximation, where in partic- 
ular the existence of oscillatory instabilities close to the linear limit can be 
proven using perturbation theory. In Section 5, SWs in the full KG model are 
analyzed numerically for two particular examples of soft (Morse) and hard 
(quartic) on-site potentials. Their linear stability properties are investigated 
using numerical Floquet analysis, and bifurcations determining the limit of 
existence for larger coupling are identified. In Section 6, we discuss the dy- 
namics resulting from the SW instabilities over short, intermediate, and long 
time scales, and Section 7 concludes the paper. The more technical parts of 
the stability analysis for the DNLS case is deferred to the Appendix. Some of 
the results obtained in this paper have been presented in an abbreviated form 
elsewhere [25]. 



2 The model 
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2.1 The Klein- Gordon system 



We consider a one-dimensional chain of classical anharmonic oscillators (Klein- 
Gordon chain) (Fig. la) with Hamiltonian 

^< + VM + -CK{un+i-u^r =E^- (1) 



n=l 



n=l 



where the total energy attributed to each site n is defined as 

= + V{Un) + ^Ck [{Un+1 - U^f + («„-! - Mn)^) • (2) 



The oscillators are submitted to a nonlinear on-site potential V{u) and are 
coupled by a harmonic inter-site coupling Ck > 0. The nonlinear potential 
can be either a soft potential such as the Morse potential (Fig. lb) 

V{u)^lu;',{e---ir (3) 



or a hard potential such as the quartic potential (Fig. Ic) 



Viu)^u;l{^ + ^). (4) 



(Recall that for a soft [hard] potential the oscillation frequency decreases [in- 
creases] with the oscillation amplitude or, equivalently, with the action.) In 
the rest of the paper we will consider V"{0) — lOq — 1. 

The choice of the Morse potential is physically justified by the fact that it 
represents an appropriate description of hydrogen bonds, and it is then espe- 
cially useful for the study of biological systems such as macromolecules and 
proteins. 

The equations of motion for the particle displacements Un then take the form 
of a discrete nonlinear KG equation, 

iln + V'{Un) - CxiUn+l + U^-l " 2'U„) = 0. (5) 



For a given solution Un{t), the corresponding linearized equations (Hill equa- 
tions) are 

e'n + ^nV"{Un) " Cx(en+1 + Cn-l " 2e„) = 0, (6) 
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Fig. 1. (a) Klein-Gordon chain of nonlinear oscillators with displacement {n}, 
coupled by a constant force Ck', (b) the soft Morse potential; (c) the hard quartic 
potential. 

where e„ denotes a small displacement with respect to the solution Un- When 
Un = (or, more generally, when V"{un) is constant), the standard harmonic 
(plane-wave) small-amplitude solutions 



e„(t) = e cos {Qn - uJi{Q)t) 



(7) 



with wave vector Q and frequency uJi{Q) are obtained. By substitution of these 
solutions in Eq. (6), we find the (phonon) dispersion relation 

LufiQ) = l + 4CKsm^Q/2 (8) 



for small-amplitude oscillations in a general potential. 



2.2 Small amplitude limit and DNLS equation 



We will now show how it is possible to describe the small-amplitude dynamics 
of general KG chains with small inter-site coupling Ck by a DNLS approxima- 
tion. Although similar results have been obtained earlier [8,16], our approach 
is slightly different and offers a deeper understanding for the role of the dif- 
ferent limits involved. Considering small-amplitude solutions to the nonlinear 
discrete KG system (5), the anharmonic terms will be small, and the on-site 
anharmonic potential V{u) can be expanded at the lowest significant order as 

V{u) = -u' + a- + f5- + ... (9) 
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For example, in the case of the Morse potential, the expansion yields a = —3/2 
and P — 7/6, while in the case of the quartic potential a = and (3 — 1. 

We then search for solutions of the nonhnear equation (5) which are locally 
close to linearized phonons (7) with wave vector Q and small amplitude e. 
If such a solution Un{t) is assumed to be time periodic with period Tf, — 
where Uh is close to some phonon frequency uJi{Q), it can be expanded as a 
Fourier series, 

Un(t)^^a^^^P^^\ (10) 
p 

where the complex coefficients fulfill = ctn~^^*- ^ timc-rcvcrsiblc solu- 
tion, the coefficients a^-* will be real. Moreover, we require the Fourier series 
to be summable as a smooth (say analytic) function of t and thus that the 
coefficients must decay asymptotically exponentially with p for each n. 

More generally, since for the linear solution the left-hand side of Eq. (5) is not 
strictly vanishing but of order e^, the perturbed solution will have a harmonic 
time dependence plus corrections of order e^. Thus, as in Refs. [8,16], we search 
for general (not necessarily time-periodic) small-amplitude solutions as series 
with the form (10) but where the Fourier coefficients depend slowly on time as 
a function of e^t. Then, in the Fourier series (10), a^^^* = a^^^^ are the leading 
order terms of order e, a^''-' and a^^-** = a^"^-* are smaller at order e^, and, since 
the Fourier components must decay exponentially, it is natural to assume that 
a^^* — a^~^) are of order e^, a^^^* = a^^^ are of order e^, etc. Then solving Eq. 
(5) at order yields ^ , with ul — 1 = 25, for each harmonic 

4°^ + 2a|a«p - CK{aT^, + ai°2, - 2af ) = (11) 
2z^,dW - 25aW + 2a(aWa(f) + a« *af ) + S/^la^paW 

-C^x(aSi + ai^li-2a«)=0 (12) 

(1 - 4c.,2)a(2) + aa« ^ - CK{a^^l + a^, 2a^) = (13) 

(1 - 9u;,2)a(f) + 2aaWa(f ) + /Ja^ ' - C^ia^ + a^^l, - 2a(f)) = 0. (14) 

The term a!^^ can be obtained exphcitly from Eq. (11) which yields 

aW = ^^^l«-H|a(i)|2 ^ith (15) 



^ Order is the lowest order taking into account the anharmonicity of the local 
potential V{u). 
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Similarly, a^^^ can be obtained explicitly from Eq. (13) but only when 2u!b > 
\/l + 4:Ck = '^z(7i"), which means that the second harmonic of the solution 
does not belong to the linear phonon band: 



,(2)=^';^^'|n-H«(l)2 ^i^h (17) 




A'= ^ (18) 

The substitution of Eqs. (15) and (17) in Eq. (12) yields a time-dependent 
extended DNLS equation with long range nonlinear interactions describing 
the slow time evolution of the dominating Fourier coefficients where the 
nonlinearities are taken into account to leading order in the amplitude e. If 
the condition 2u!b > y/l + ^Ck is not fulfilled, \rj'\ = 1 [j]' complex, on the 
unit circle) and the sum (17) defining the interaction range in general diverges. 
Thus, for a general potential with non-zero a, this condition imposes an upper 
validity limit for the coupling strength Ck for this approach. (In the linear 
limit ujb = uJi{Q), the condition is fulfilled for all Q only when Ck < 3/4.) 

In the particular case when a = (e.g. the quartic potential), there are no 
long-range terms in Eq. (12), Q which immediately reduces to the ordinary 
DNLS equation with only nearest-neighbour interactions (Eq. (19) below). 
For a general potential, we consider the case of small Ck (i-e. close to the 
anticontinuous limit Ck = 0). Then, — 1 = 26 will be small of order of Ck, 
and v = Ck + 0{Cl), r]' = -^ti + 0{Cl). Considering the modified DNLS 
equation (12) at order Ck and neglecting terms of order e^Ci^, the long-range 
terms disappear and it reduces to the standard DNLS equation which well 
describes the dynamics of small-amplitude oscillations in KG chains close to 
the anti-continuous limit: 

ilpn = Sljjn - a\lpn\^i)n + C(^n+1 + " ^ijjn), (19) 



where = y^a^, A = -5a'^/3 + 3/3/2, a = sign(A) and C = Ck/2. 
Substituting the corresponding values of the parameters a and (3, we obtain 
that for soft potentials A < and then a = —1, while for hard potentials A > 

^ When a = 0, coupling to higher harmonics appears only to order and higher. 
In particular, for the quartic potential only odd harmonics couple to each other. 
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and a — 1. Actually both cases cr = ±1 arc formally equivalent at this order 
of the expansion, since by changing i/jn into (— I)"-?/'* we get a new equation 
for {ipn} which has the same form as Eq. (19) but with S replaced by (4C — 5), 
a by —a and C is unchanged. So within the DNLS approximation we can 
assume for convenience a — —1, corresponding to a soft potential, without 
loss of generahty. 



3 Standing waves 

3.1 Standing waves represented as map trajectories 

SWs in the nonlinear KG model (5) can be searched as time-reversible and 
time-periodic solutions {it„(t)} with frequency oub as in Eq. (10). The Fourier 
coefficients w^^"' of the time-periodic function V'{u{t)) = w/^^ e'*"^''* define 
nonlinear functions {w^^\{a'^'^^})} of the Fourier coefficients {a*^*^} of u{t). 
Then, Eq. (5) can be written as 

- + ^^W}) - CKia^S + «Si - 2a^)) = 0. (20) 

This equation defines a symplectic map 

iWllih W^n) = '5({a(f)}, {a^fl,}), (21) 

where the SWs should appear as bounded trajectories. Since this map operates 
in an infinite dimensional space, its study is quite uneasy and moreover the 
numerous theorems about symplectic maps available for finite systems do not 
hold for infinite systems. It is convenient to reduce the dimensionality of this 
map, which can be done using the DNLS approximation derived in Sec. 2.2 
(generally valid in the limit of small- amplitude solutions and small inter-site 
couphng) to ehminate all with p 7^ 1. 

Since the SWs are time-periodic solutions, all the amphtudes a^^ of the Fourier 
expansion (10), and therefore all the ipn, will be time independent, and they 
can be chosen real by restricting to time-reversible solutions. Then, vanishing 
the left-hand side of Eq. (19) yields 

Si^n - (J^l^l + C(^„+i + - 2V^„) = 0. (22) 

This equation defines a nonlinear symplectic cubic map S in the 2-dimensional 
real plane, which recursively determines {il^n+iiil^n) from {il^m'^^n-i)- With the 
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change of scale t/jn — VCi/^'n, this map depends only on the parameter S' = 
S/C: 

ij'^^, = (2 - + ai^'^ - , < = (23) 



As is well-known [24], this map exhibits a rich variety of orbits, including 
elliptic and hyperbolic fixpoints and periodic cycles, KAM tori, Aubry-Mather 
Cantor sets (cantori) [15,13,26] and chaotic orbits (see Fig. 2). Since we are 
searching for nonhnear SWs in continuation of the linear SWs, only bounded 
trajectories which can be continued to zero amplitude by varying S' — are 
of interest. 
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Fig. 2. The symplectic map S (oi S' = ^ = 0.5 (cr = —1). The visible chain of 
alternating elliptic and hyperbolic periodic points correspond to the two families of 
SWs with Q = 7r/3. 

When 6 belongs to the phonon band, that is when < 6 < 4C, the lineariza- 
tion of Eq. (22) has the solution ipn = cos{Qn + (p), where Q fulfills 



5 = 4Csin2Q/2 = 5o(Q), 



(24) 



which is equivalent to the KG dispersion relation o;^ = ujf{Q) ~ l-|-4Cft:sin^ ^ 
(Eq. (8)). Then, the trivial fixpoint Fq = (0, 0) of the map S is hnearly stable, 
that is elliptic with a rotation angle Q. This solution corresponds to a linear 
SW and the rotation angle is nothing but the wave vector Q of the SW. When 
6 > Sq{7i) = 4C the fixpoint Fq = (0,0) becomes hyperbolic with refiection, 
while it becomes hyperbohc without refiection for S < So{0) — 0. 
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When (V = > = 4 and a = —1, it can be proven that the fixpoint 
Fq = (0, 0) corresponds to the only bounded trajectory of the map and that 
all other trajectories {ipn} go to infinity for either positive or negative n or 
for both. To demonstrate this assertion, let us assume that at some point a 
trajectory fulfills the condition 7^ > Then Eq. (23) yields that 

\^Ul\>i\'^-^'-€\-mnl (25) 



which for S' > 4 yields > (5' — 3)|'(/''J > [■(/''J and thus implies recursively 

that {ipnl goes to infinity for large positive n. If we have the reverse inequality 
lipnl ^ I'V'n-il; wc obtaiu identically the same result but for n going to —00. 
Consequently, there are no spatially bounded time-periodic solutions with a 
frequency above the phonon band for a DNLS equation with a soft potential 
{a = —1). For a hard potential {a = 1), the same argument excludes non- 
trivial bounded time-periodic solutions below the phonon band 6' < 0. 

When 5' = |; < 4 for soft potentials (or, equivalently, 5' > for hard poten- 
tials), it can also be proven that the trajectories which escape too far from 
the origin go to infinity. Indeed, for a = —1, if ip'^ > 4 — S' and lipnl > |V'n-i|, 
the inequality (25) implies that ipn+i > "^n > 4 — 5' and recursively that \il)n\ 
goes to +CXD for n — > +00. A similar result holds for the reverse inequality 
\'4^n\ ^ |V^n-i| but then j-j/'^l goes to -|-oo for n — > —00. Nevertheless, when 
5' < 4 the nonlinear map S exhibits many bounded trajectories representing 
SWs. As mentioned above (see Fig. 2), some of them are chaotic but there 
are also periodic cycles, KAM tori and cantori for example around the elliptic 
fixpoint Fo = (0,0). 

As the nonhnear SWs we search for should be continuable to zero amphtude 
when varying 5' = ^, the corresponding orbits should merge into the elliptic 
fixpoint -Fo for the particular value 5 = So{Q)- As suggested in [12], trajectories 
which have this property should be either periodic cycles or quasiperiodic 
trajectories. The quasiperiodic trajectories, corresponding to wave vectors Q 
incommensurate with 27r, should, for most irrational values of Q/^tt, appear 
as KAM tori with rotation angle Q. Thus, given a rotation angle Q, there is 
at least one trajectory which appears with this rotation angle at the elliptic 
point Fo for S = So{Q), and which for 6 close to 6o{Q) is located close to Fq. 
This trajectory, which for rational Q/27T = r/s is a periodic cycle with period 
s and for irrational Q/Stt generally a KAM torus, represents a small amphtude 
SW with wave vector Q. For a general potential (soft or hard), the coefficient 
a of the cubic term in the DNLS equation (22) determines the behaviour 
of the tori around the elliptic fixpoint Fq. For a soft potential {a < 0) the 
rotation angle for orbits rotating around Fq increases with increasing radius, 
and in particular for KAM tori with small radius the rotation angle can be 
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calculated to lowest order in a^ 



cosQ^l-| + ^. (26) 



Since the map S depends continuously on the map parameter S', the periodic 
cycles, KAM tori and cantori of S can be continued as a function of 6', and its 
orbits at fixed rotation angle depend continuously on 6' [12]. Then, as So{Q) is 
a monotonously increasing function of Q, we can conclude that SWs with wave 
number Q exist only for S < So{Q). Similarly, for a hard potential (a > 0) 
the rotation angle decreases with increasing radius and SWs exist only for 
S > 6oiQ). 



Q/ 2jr rational (Commensurable SWs) 




(H) hyperbolic without reflection 




(E) hyperbolic . , -. elliptic ^-^^^^^ 
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Q/2% iiTational (Incommensurable SWs): 




(H) Cantorus 




^^^^^ KAM torus * 




(E) midgap 
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8 5 





Fig. 3. Evolution of the trajectories for variable b' and fixed wave vector Q [a = — 1). 



To study the evolution of the trajectories representing SWs (for a = —1) 
when varying the map parameter S', we thus start from the linear limit at 
S = So{Q) corresponding to the fixpoint Fq and decrease the value of 6. This 
evolution is schematically illustrated in Fig. 3. For a rational value of ^ = r/s 
(r and s irreducible integers), which corresponds to SWs commensurable with 
the lattice periodicity, there are two families of SWs corresponding to the two 
s-periodic cycles which bifurcate in pair from Fq at the linear limit 6 = Sq{Q). 
We will call these s-periodic cycles 'h-cycle' and 'e-cycle', respectively. The 
/i-cycles are hyperbolic without reflection for all 6 < 6o{Q). By contrast the 
e-cycles are elliptic trajectories for 6 close to 6o{Q), but decreasing the value 
of 6 they become hyperbolic with reflection below a critical value 6 = Sc{Q) 
which can be calculated accurately numerically e.g. using Greene's residue 
[27]. 

When Q/2tt is a generic irrational number (i.e., not a Liouville number), the 
trajectories representing the spatially quasiperiodic, incommensurate, SWs 
emerge from Fq as the KAM torus with rotation angle Q, which is unique 
apart from a phase shift. Decreasing the value of 5 below the critical value 

^ This follows e.g. from Eq. (A. 13) in the Appendix using Eq. (43). 
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Fig. 4. (a) and (b): Hull functions xs{x) for SWs in a DNLS model (22) with 
Q/27T = 233/610 ~ (3 - V5)/2 = ac and (a) d' = 3.0 resp. (b) 6' = 2.8 (C = 1). 
The analyticity is found to be broken at 6' ~ 2.92 for Q/2Tr = ctg- (c) and (d): 
Hull functions Xk{x) for SWs in a KG chain (5) with a Morse potential, with 
Ck = 0.05, Q/27r = 34/89 ~ gg and (c) uj^ = 1-072 resp. (d) uj^, = 1.065. The 
breaking of analyticity occurs at Wf, ~ 1.070 for Q/2'i: = ac- 

S = Sc{Q), the KAM torus bifurcates into two types of trajectories, the can- 
torus and its associated 'midgap' trajectory [28]. This midgap trajectory con- 
sists of isolated, nonrecurrent points in the gaps of the cantorus with rotation 
angle Q, to which it is homoclinic. As the KAM torus is an analytic trajectory 
while the cantorus and the midgap trajectory are non-analytic, the transition 
at the critical point 6c corresponds to a transition by breaking of analyticity 
(TBA) [13,26]. The transition point 6c{Q) is determined as the limit of the crit- 
ical values 5c (Q*-^^) for a sequence of periodic e-cycles whose rotation angles 
Q^'^'* are obtained from successive rational approximants of Q/2-k [27]. With 
this transition the two families of non-analytic incommensurate SWs, repre- 
sented by the cantori and the midgap trajectories, respectively, merge into the 
unique family of analytic 5*1^ defined by the KAM tori for 5 > SdQ). Defining 
the hull (or envelope) function Xsi^) for the SW as 

^n = Xs{Qn + (P), (27) 

which gives the shape of the wave rescaled inside one period 27r, the TBA for 
this hull function when decreasing 6 is shown in Figs. 4(a) and 4(b). As the 
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model symmetry implies two axis of symmetry for the KAM tori (see Fig. 2), 
the hull function Xsi^) can be chosen with the symmetries Xsi^) = ~Xs{~x) 
and Xs{x) = ~Xs{x + tJ")- The nonanalyticity of the SW hull functions for 
large negative 6 makes a description in terms of multibreathers natural, as will 
be discussed in Sec. 3.2. 

As the trajectories of the map S represent approximate solutions for the KG 
model in the limit of small coupling, we should expect to find this TBA also 
for SWs in the KG system (5). We can calculate numerically exact KG SW so- 
lutions using a Newton scheme as described, e.g., in [11,7,29] (see also Sees. 3.2 
and 5). Defining the hull function Xk{x) for these SWs as ■u„(0) = XK{Qn + (p) 
and continuing these solutions versus Ub to zero amplitude, a TBA is found for 
the function Xxi^) as illustrated in Figs. 4(c) and 4(d). (Note however that for 
a non-symmetric potential such as the Morse potential, the hull function does 
not in general have the above-mentioned symmetries, i.e., Xk{x) 7^ ~XKi~x) 
and XKix) ^ -Xk{,x + tt).) 



3.2 Standing waves as multibreather solutions 

When the SW frequency LOb is below the phonon spectrum for a soft potential 
(that is ujb < 1) or above for a hard potential {ujb > y/l + 4:Ck) and close 
enough to the anticontinuous limit Ck = 0, the SW solutions can be iden- 
tified as multibreathers [4] with period Tb = These are oscillatory states 
composed by nonlinear combinations of severalj^ time-periodic and spatially 
localised breathers in a single nonlinear system [4,30]. These multibreather 
solutions are characterized by a coding sequence a„ [4,12] which defines the 
particular oscillatory state of the oscillator n in the case of zero coupling. As 
for general time-reversible time-periodic solutions with period Tb each indi- 
vidual oscillator can oscillate with frequency joob {j integer) with two possible 
choices for the phase, we can, at Ck = 0, characterize all such solutions for 
the whole system by the following coding sequence: 

• 5"„ = when the oscillator n is at rest; 

• 5"„ = +1 when the oscillator n oscillates with a finite amplitude at frequency 
Ub and in phase with the chosen reference solution of the single oscillator; 

• 5"„ = — 1 when the oscillator n oscillates at frequency Ub in antiphase with 
the chosen reference solution; 

• an = +2 when the oscillator n oscillates in phase at frequency 2u!b] 



^ We will use the term multibreather to refer to combinations of either finitely or 
infinitely many individual discrete breathers. 
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Fig. 5. Continuation of the multibreather solution with wave vector Q into the hnear 
phonon band for a soft potential, from the anti-continuous limit Ck = to the linear 

limit ijJh = y''i + '^Ck sin^ ^ illustrated by the gray line. (For a hard potential, the 
continuation starts at Wft > 1 and can be performed at constant frequency.) 

These solutions are proven to be continuable when Ck varies from the an- 
ticontinuous limit up to a finite value of the coupling. At this value, which 
depends in general on the solution in consideration, bifurcations are found and 
the solution can be lost. This will be numerically investigated in more details 
for SWs in Section 5. We consider now only the continuable regime of the 
multibreathers, and in particular we search among these multibreather states 
the solutions which are in continuation of the linearized SWs at low amplitude. 
For that purpose, we have to find an appropriate coding sequence such 
that the multibreather which exists at Ck 7^ can be continued versus its fre- 
quency inside the phonon band, as illustrated in Fig. 5. A similar problem has 
already been studied in [11] for a disordered model. Actually, it was found that 
most of the multibreather solutions undergo bifurcations and disappear before 
their frequency enters the phonon spectrum, although some of them can be 
continued while their frequency belongs to the phonon spectrum. For this last 
type of solutions an ansatz was found for calculating specifically the multi- 
breather states which are continuable as linearized harmonic modes. However, 
the method used in [11] was explicitly based on the property that the phonon 
spectrum is discrete with localized eigenmodes. In our translationally invariant 
model this property is not fulfilled so that we cannot use the ansatz proposed 
in [11]. There are nevertheless two linearized solutions which can be trivially 
obtained by continuation of multibreather states, namely the Q = phonon 
from the multibreather with code ... + 1 + 1 + 1 + 1 + 1 + 1 + 1 + 1..., and the Q = vr 
phonon from the multibreather with code ... + 1 — 1 + 1 — 1 + 1 — 1 + 1 — 1.... Our 
problem is to find the ansatz which gives the coding sequence for a general SW 
with wave vector Q. Once this code is obtained, it is possible to calculate the 
nonlinear SW for different values of the frequency and coupling by numerical 
continuation. 

The method we propose to find the ansatz consists in passing from the KG 
model (5) to the DNLS model (19) (where we for simplicity assume a = —\ 
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without loss of generality) by using the small-amplitude (and generally also 
small-coupling) approximation (Sec. 2.2), so to obtain the SW representations 
as trajectories of the symplectic map S (23). As the DNLS map S only depends 
on the parameter S' = the anticontinuous limit C = for fixed frequency 

5 < is in the DNLS description equivalent to the large-amplitude limit 

6 — oo for fixed coupling C. Thus, continuing the SW trajectories to this 
limit, they can be described as multibreathers and associated to a coding 
sequence. As for C = the stationary solutions to Eq. (22) can only take the 
values ipn = ±a/— 5 and ipn = 0, the code will be defined as dn = +1 when 
ipn = \/—^, O'n = —1 when ip^ = and = when ipn = 0. These 
codes can then be used as initial conditions in an iterative Newton scheme to 
calculate numerically the SWs in the original KG model (5) to any desired 
accuracy. 

By continuing the SW solutions calculated for the DNLS map to large negative 
6' we guessed and tested, by inspecting the map trajectories, the following 
ansatz for their coding sequencesj^ 

&n = Xo{Qn + (j)), (28) 



where o"„ can be quasi-periodic or periodic, depending on whether Q/2n (with 
< Q < it) is irrational or rational, and the hull function Xo{x) is a 27r- 
periodic odd function defined for x G [— tt, it] as 



1 for (tt - Q)/2 < a; < (tt + Q)/2 

-1 for - (7r + g)/2 < X < -(tt -Q)/2 

elsewhere 



(29) 



This function Xoi^)^ rescaled by a factor y/—6, is the limit for C = of the 
hull function Xsi^) introduced in Section 3.1. 

This method for generating the code can be visualised by a graphical circle 
construction (see Fig. 6). To construct the code for a nonlinear SW with wave 
vector Q, the circle is divided in four parts to which the values 0, +1 or —1 are 
associated. Rotating around the circle by the angle Q yields the code sequence 
of the SW. The value of the phase 0, corresponding in the circle construction 
to the initial point for the rotation, has particular importance. For all phases 
(f) ^ (j)^ = ±(7r — Q)/2 — mQ {m integer), the coding sequence generated does 
not contain any consecutive codes -|-1 or —1. We call these SWs 'type H' and, 
with the terminology used in Sec. 3.1, they correspond either to the /i-cycles 

^ For a hard potential {a = +1), Eqs. (28)-(29) yield the SWs with wave vector 
TT — Q via the transformation dn (— l)"'^^. 
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Fig. 6. Graphical representation of the generation of the code sequence for a general 
nonlinear SW with wave vector Q. 

(for rational Q/2ti) or to the cantori (for irrational Q/2tx). For the particular 
phases = 0^, x — Qm + is at a discontinuity point of Xq{^)-i ^^^d, as 
a consequence, cr„ has two consecutive +1 (or —1). These SWs, called 'type 
E', correspond to the e-cycles (for rational Q/2ti) or to midgap trajectories 
(for irrational Q/2t{). As will be seen in Sec. 4, this characteristic of the code 
will imply important consequences for the stability of nonlinear SWs. We note 
that the coding sequences for SWs of type E will always be symmetric around 
a bond center, and can therefore be chosen to fulfill cr„ = cr_„+i by choosing 
(f) = {'K — Q)/2. In the special cases when Q = 2k'+i ^ ^' integers), the codes 
for type E SWs will in addition be antisymmetric around some lattice site, and 
can therefore alternatively be chosen to fulfill cr„ = — (J-n by choosing (p — 0. 
The coding sequences for type H SWs are antisymmetric around a lattice site 
except when Q — Ifc'+i ^' which case they are antisymmetric around a bond 
center. In the former case we can choose = so that a"„ = —cr-m while in 
the latter case 5"^ = —a_n+i by choosing = —Q/2. 

Finally, we remark that the construction (29) to generate the SW coding 
sequence can, for type H SWs with Q close to vr, be regarded as a periodic or 
quasiperiodic repetition of defects or discommensurations added to the Q = n 
phonon, where each discommensuration consists of inserting one additional 
site with code ct^ = in the coding sequence for the Q — ir phonon. In the limit 
Q — > TT, the distance between the discommensuration sites goes to infinity, so 
that we recover the case of one single defect site which was studied for the 
DNLS model in [18] and shown to correspond to the discrete counterpart of 
the dark NLS soliton. As will be discussed later, this point of view is useful 
in order to understand the origin of the instabihties of the SWs for Q > 7r/2. 



17 



4 Stability of nonlinear standing waves in the DNLS approxima- 
tion 

We now turn to the main topic of this paper, namely the characterization of 
the dynamical stability properties of nonlinear SWs. We will here first consider 
small-amplitude solutions to the KG equation (5) with small Ck, so that the 
DNLS approximation (19) is well justified. In the next section we will study 
numerically the original KG model for larger coupling to investigate the limit 
of validity for the DNLS approximation. Here, we assume for convenience a = 
— 1 without loss of generality (Sec. 2.2). With the substitution ipn 'ipn+^nit), 
where e„ is small, the linearization of Eq. (19) yields a standard Hill equation 
for en [31], 

iin = Sen + '2,\i)n\^en + i^lel + C(e„+i + e„_i - 2e„). (30) 

A solution {V'n} of Eq. (19) is said to be linearly stable when there is no 
solution {e„(t)} of Eq. (30) which diverges exponentially at large time. 

In the case when ipn is a real and time-independent solution of Eq. (19), the 
general solution of (30) can be searched as linear combinations of solutions 
with the form 

en(t) = a„e-*+6;e--* (31) 

where the eigenfrequencies uj = ujs and eigenmodes {a„, 6„} are determined 
by the eigenequation 

(2C -5 - 2ipl)an - C(a„+i + a„_i) - iplb^ = Us^n 

il^lan + C{hn+i + 6n-i) - (2C - 5 - 2,pl)hn = ^.6„. (32) 

Then, we note that if Ug is an eigenvalue for the eigenvector {a^, 6n}, then —Ug 
and cj* are eigenvalues for the eigenvectors a„} and {a* , 6* }, respectively^. 
As a result, linear stability is achieved if and only if all the eigenvalues of Eq. 
(32) are real. Thus, for a large but finite system with sites, we should have 
2A^ real eigenvalues ujg for a linearly stable solution. 

It is convenient to redefine new variables = a„ + 6„ and Wn = an — bn, so 
that Eqs. (32) become [31] 

LoWn = {2C-6- iPDWn - C{Wn+l + Wn-l)=0OgUn 
C,Un = {2C-6- SijDUn - C{Un+l + Un-l)=UjMn. (33) 

^ Tlie exponential of this operator is symplectic. 
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The operators Co and Ci are self-adjoint, but not the whole operator C = 
Co\ 

. Note that cu^ = is always an eigenvalue of C for f/„ = 0, since the 

A o) 

first equation in (33) then has the explicit solution Wn = ipn, corresponding 
to the time-independent phase mode e„ = iipn resulting from the invariance of 
the model under global phase rotations. Moreover, for an incommensurate SW 
with analytic hull function (27) there is also an additional solution to (33) with 
Us = and Wn = 0, namely the sliding mode solution f/„ = dxs/d(l> to the 
second equation, resulting from the invariance of analytic SWs under spatial 
translations. See also Appendix A. 2 for a further discussion of the spectrum 
of the operators Cq and Ci. 

We also define the Krein signature [32,4] associated with a pair of eigenvalues 
icus as [33] 

K,{uJs) = sign^ [|aj2 _ |^^|2] ^ sig^ ^ [f/^iy^* + w^K] , (34) 



where {a^, 6„} and {[/„, Wn} are the eigenvectors of (32) resp. (33) with eigen- 
value +uJs- In general, when varying the parameter 5', instabilities may occur 
only through collisions between eigenvalues with different Krein signatures. 

When 5 < and C — > (anticontinuous limit), ipn = — 5 or (codes dn = ±1 
or 0, respectively), and all eigenfrequencies of (33) are easily seen to be real. 
For a finite system with P sites where |5"„| = 1 and Z = N — P remaining 
sites where dn = 0, there are 2P eigenvalues degenerated at Ug = 0, while 
the remaining eigenvalues are located as Z pairs at a;^ = ±5. Increasing the 
value of the coupling, the eigenvalues will move losing their degeneracy and, if 
a pair of eigenvalues goes out in the imaginary plane, the nonlinear SW may 
become unstable. In particular, when C is increased from zero, the eigenvalues 
at ujg = ±6 (which all have Krein signature +1) will move only along the real 
axis causing no instabilities for the SWs. On the other hand, the behaviour of 
the eigenvalues originating from a;^ = will be different for the E'-type and H- 
type SWs (as defined in Sec. 3.2). As the coding sequences for the E-tjpe SWs 
contain pairs of consecutive +1 (or —1) and the anharmonicity is soft, it follows 
from a general result [4] that some pairs of eigenvalues originally at 0;^ = 
leave the real axis and move out on the imaginary axis as soon as C 7^ 0. 
As a consequence, the corresponding SWs become dynamically unstable for 
6 < and small C. Moreover, we find numerically that this instability persists 
for all 6 < So{Q), so that SWs of type E are unstable for any nonvanishing 
amplitude.^ By contrast, for H-tjpe SWs, for every site with code 5"„ = 1 the 

^ This non-oscillatory instability is also proven in the Appendix for commensurate 
SWs of type E in the small-amplitude limit. 
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nearest surrounding sites with nonzero code have cr„ = —1, and reciprocaUy 
(i.e., there is a phase shift of tt between any two consecutive sites with nonzero 
oscillation at C = 0, neighboring or not). Then, the theory of effective action 
developed in Ref. [4] yields, with similar calculations as in [34,18], that all 
eigenvalues originating from cus — (having Krein signature —1 opposite to 
those originating from oug — ±5) move away from the origin along the real axis, 
except one pair remaining at cjg = (with Krein signature 0) corresponding 
to the global phase invariancc. Thus, the stability of these SWs is preserved 
for C not too large (c.f. Fig. 7 (a)). 



4.1 Commensurate SWs 



For commensurate SWs with Q = nr'/ s' (r' and s' irreducible integers) in an 
infinite system, the lattice periodicity of ip'^ in (33) implies that the degenerate 
eigenvalues on the real axis form bands as C is increased from zero (see Fig. 7) . 
As there are s' — r' sites with code o-^ = in a lattice period of s' sites, there 
will be s' — r' symmetric pairs of bands with Krein signature +1 originating 
from cus — and separated by gaps. The remaining eigenvalues form r' 
symmetric pairs of bands with Krein signature —1 originating from Ug = 0. 
As the eigenvalue a;^ = always remains in the spectrum, there is no gap at 
zero which thus always remains at a band edge. For fixed C and increasing S, 
the bands originating from ci;^ = and those from cug — ±5 broaden and the 
gap separating these two classes of bands will shrink. At a critical value S ~ 
^k{Q) this middle gap shrinks to zero and, as shown in Fig. 7(b), symmetric 
pairs of eigenfrequencies with opposite Krein signature start to collide and 
move out from the real axis into the complex plane. Thus, oscillatory (Krein) 
instabilities (i.e., caused by complex eigenvalues) of the nonlinear SW are 
induced. 
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Fig. 7. Eigenvalues of Eq. (33) for a type H SW with Q = 37r/4 and (a) 5' 
(b) S' = -7.5, (c) S' = -4.5 resp. (d) S' = -0.5 (C = 1). 
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This instability mechanism is similar to that previously observed [18] for the 
'discrete dark solitons', where the anticontinuous coding sequence consists of 
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one single zero-amplitude 'defect' site added to a background wave with Q = ir 
(the latter being always stable [8]). In fact, as was mentioned in Sec. 3.2, for 
Q close to TT the SWs generated from (29) can be viewed as a periodic (or 
quasi-periodic for incommensurate SWs) repetition of defects or discommen- 
surations added to the stable wave with Q — n. Thus, the SW can be regarded 
as a lattice of weakly interacting discommcnsurations ( "dark solitons" ) , each 
of which shows an oscillatory instability due to a resonance between a trans- 
lational mode localized at the discommensuration and extended modes [18]. 
This oscillatory instabihty detected numerically for the DNLS model can also 
be confirmed by a direct numerical Floquet analysis of the KG Eq. (5) (see 
Section 5). 

Increasing further 6 above the critical value 6k {Q), bands continue to overlap 
and generate new instabilities. For many commensurate values of Q, it is found 
that when S increases beyond Sk{Q), the SW always remains unstable until it 
vanishes with zero amplitude at the linear limit d = ()q{Q). However, for some 
rational values Q/27r = - with small s, in some intervals of 6 the band overlap 
disappears, and consequently the stability is temporarily regained (e.g. for 
Q = 37r/4, the H-type SW is stable for -4.9 < S' < -3.9, see Fig. 7(c)). In 
any case, for all commensurate SWs with Q ^ tt, there will always be a final 
interval Si{Q) < S < Sq{Q) where bands with opposite Krcin signature overlap 
(Fig. 7). As a consequence all commensurate SWs with Q ^ n are unstable 
for small but nonzero amplitude. This result can be simply understood from 
the remark that at the linear limit 5 — Sq{Q) where ip'^ — 0, the frequencies 
of the (gapless) spectrum of (33) are given by 

Go's = ±2C| cos Q — cos g|, (35) 

< g < 271, and the corresponding Krein signatures are given by sign(cos(5 — 
cosg). Thus, in the linear limit there is an interval around Us = where 
eigenvalues with opposite Krein signature overlap, although there is of course 
no instability since there is no coupling between modes with opposite Krein 
signature in this limit. By a continuity argument, since in the commensurate 
case when ip"^ ^ a finite number of bands with nonzero widths appear for 
5 < So{Q), bands with opposite Krein signature would, assuming the absence 
of mode coupling, continue to overlap on the real axis close to ci;^ = until S 
reaches some value 5i{Q) which must be strictly smaller than So{Q). However, 
as in general there is coupling (i.e., resonances) between modes of different 
Krein signature, the overlapping bands move out in the complex plane and 
instabilities result. This is proven in the Appendix using perturbation theory 
and the method of 'band analysis' outhned in Sec. 4.2 below. 

Let us also remark that in the case of finite-size systems (e.g., with periodic 
boundary conditions), the eigenvalues ujs will evidently be separated in a dis- 
crete spectrum without forming continuous bands. The initial scenario when 
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increasing 5' from — oo will however be qualitatively the same as described 
above, and oscillatory instabilities will result similarly as illustrated in Fig. 
7(b) (although the bands are replaced by discrete eigenvalues). However, the 
above argument showing instability close to the linear limit for infinite sys- 
tems is obviously not valid for finite systems, as there are no continuous bands 
around Ug = Q but only discrete eigenvalues. Thus, overlap (resonances) can be 
avoided and, in fact, we find that for finite systems stability will generally be 
regained in some interval close to the linear limit. The scenario by which the 
stability is recovered is through 'reentrant instabilities' similarly as in [17,18], 
and appears in an analogous way also for the KG model as will be described 
in Sec. 5. 

4.2 Incommensurate SWs 

For incommensurate SWs [Qjl-n irrational), oscillatory instabilities will occur 
at some critical value ^Kiff) just as for the commensurate case described 
above, although the spectrum of (33) now has a Cantor set structure with 
infinitely many gaps. However, for small-amplitude (analytic) incommensurate 
SWs, the Cantor set nature of the spectrum (with 'bands' of zero width) 
implies that the argument used in the previous section cannot be used to prove 
the instability of the nonlinear SWs close to the linear limit, since it would only 
yield 5i(Q) = 5o(Q)- Indeed, numerical analysis of finite-size approximants to 
incommensurate analytic SWs with \Q\ > n j1 typically yields a scenario with 
reentrant instabilities which become very weak and practically invisible for b 
close to 5o(Q)! and extrapolating this scenario to infinite systems could lead 
to the (erroneous) guess that analytic SWs with \Q\ > 7i/2 could be stable for 
some 6 (possibly constituting a Cantor set as observed in disordered systems 
[11]). This scenario would mean that in the absence of mode coupling, the 
(Cantor-like) subspectrum of (33) with positive Krein signature would, for 
some 6, exactly fit into the gaps of the other (Cantor-like) subspectrum with 
negative Krein signature, and thus instabilities could be avoided. However, 
as 5 — > 6o{Q) both subspectra acquire full measure, and if their measure is 
assumed to be continuous with respect to 6 overlap cannot be avoided for 
6 close enough to 6o{Q), and instabilities should result also in this case. In 
fact, this intuitive argument is confirmed by the more sophisticated analysis 
described below (some technical details are carried out in the Appendix), 
which shows that also the incommensurate SWs are unstable for 6 close to 
So{Q), but when \Q\ > 7r/2 only through higher order instabilities, which are 
thus more difficult to detect numerically or experimentally. 

To prove the instability of analytic SWs we use here the method of 'band 
analysis' [4] .[^ Its basic idea is to embed the non-Hermitian eigenvalue problem 

^ The same approach can be used for small-amplitude commensurate SWs, see 
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(33) into a wider eigenvalue problem but for a Hermitian operator, which has 
more properties. Then, according to [4], we associate this problem with an 
extended eigenvalue problem 

ien - Sen - 2|^/'npen - i'l^-n " C{en+i + e„_i - 2e„) = Ee„. (36) 



The eigenvalues and the corresponding eigensolutions of this eigenequation can 
be searched again with the form (31), where a; is a real parameter. Then, the 
eigenvalues E with eigenvectors {a„, 6„} are those of a self-adjoint operator: 



{-u + 2C - 5 - 2i)l)an - C(a„+i + a„_i) - = Ban 

-i^lan - C(6„+i + hn-i) + {uj + 2C-5- 2^^)^^ = Ehn- (37) 

Writing as before = a„ + 6„ and Wn = an — bn, (37) becomes 

{2C-6- 3^l)Un - C{Un+l + Un^l) - OoWn = EUn 

-uUn + (2C - 5 - i;l)Wn - C{Wn+l + Wn-l) = EWn, (38) 



or 

-00 Co I \ {Wn} I \ {Wn} 



A i{u„) ^^^^ 



with Cq and Ci defined by Eqs. (33). Then, each eigenvalue Ei,{u) of the 
problem (39) becomes a smooth, real function of u called band (symmetric 
with respect to = 0), which also depends smoothly on the other model 
parameters (as before, we only need to consider variations of the parameter 
5' = ■^). If we assume that the parameter ui represents the eigenfrequencies 
for the initial non-Hermitian problem (33), the problem (39) reduces to the 
problem (33) for the particular value E = 0. Thus the real eigenfrequencies 
Us of (33) are determined by the intersections E^ius) = of the bands with 
the axis = 0, and their Krein signature by sign(— ^^) aX u = Ug > Q [4]. 
Varying the model parameter 5' the bands will generally move in the (cj, E)- 
plane, and if a band loses a pair of intersections with the zero-energy axis (i.e., 
if a 'gap' opens around this axis) an instability occurs for the initial problem 
(33). □ 



Appendix. 

^ For an oscillatory (Krein) instability, two symmetric pairs of intersection are 
simultaneously lost. 
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At the linear limit {ijjn — > and 5 — > So{Q)), the eigenvalue problem (39) is 
easily solved by introducing a wave vector q and writing 

U{q) =T.Un e'^^ W{q) = E e'^" . (40) 

n n 

At this limit, there are two sets of plane wave solutions [/„ = ±H^n — e'^" with 
corresponding eigenvalues 

Eo±{q; a;) = 2C - So{Q) - 2Ccos g ^ a; = 2C(cos Q - cos q) ^ lo. (41) 

For a finite system with N sites, = 27rk/N with < k < N integer, and 
thus there are 2A'^ bands described by straight lines with slope =f1 and indexed 
by the wave vector q (Fig. 8 (a)). Each band Eo±{q;uj) intersects E = at 
u!±{q) = ±(2C — do{Q) — 2C cosg) = Ws{q), so that we recover the result (35) 
that Eq. (33) has 2A'" eigenvalues on the real axis so that the solution ■^,1 = 
is linearly stable as expected, and remains stable also in the limit of large N. 

When ipn is small but non-zero, it can be treated as a perturbation for the 
operator in the left-hand side of Eq. (39). Then, gaps will generally open in 
the continuous band spectrum of this operator at points or lines of degeneracy. 
Such a gap might suppress a band intersection Eq±{u!) = with the real axis, 
and thus generate a linear instability. There are two kinds of degeneracies. The 
first kind occurs for both branches of eigenenergies EQ^{q; u) and EQ_{q; uj) for 
all UJ, since Eo±{q; uj) = Eq±{271 — q; uj)). The second kind occurs at transverse 
intersection points, where bands with opposite slope intersect if Eo_{q;uj) ~ 
Eo+{q';u!) for some q, q'. 

For an analytic SW, the hull function (c.f. Eq. (27)) Xsi^) V'n tie chosen 
even and 7r-pcriodic, and thus can be expanded in a Fourier series as 

V'^=E/pe''"^"'^^^^ (42) 

p 

where fp = f_p. We can choose the arbitrary phase as = 0. Since, as for 
conventional perturbation theory, the total sum of the gap openings will be 
proportional to Y,p \fp\, perturbation theory will be consistent in principle only 
when this sum is convergent and small. Thus, perturbation theory is applicable 
for incommensurate SWs only when Xsi^) is analytic, in which case fp decays 
exponentially for p — > oo, and not when ip^ is represented by cantori because 
then this sum is divergent. Then, considering only the regime close to the 
linear limit where the SW is analytic with small amplitude e^, Xs{x) becomes 
close to a harmonic function cos x, and we have 

/o ^ 2/i = 2/_i ^ e1/2 > 0, (43) 
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while the other Fourier coefficients fp will be small of order e^l^L Using the 
space Fourier transform (40), Eq. (38) becomes without approximations 



(2C-5- 2Ccos q - 3/o)C/(g) - 3 fpU{q + 2pQ) - uWiq) = EU{q) 
-uUiq) + {2C-6-2C cos q - fo)W{q) - ^ fpW{q + 2pQ) - EW{q). 

(44) 

Thus, this equation couples the space of solutions {U{q),W{q)) with wave 
vector q of the unperturbed problem (39) with those with wave vector q + 2pQ, 
and the coupling strength is proportional to fp. The matrix of the operator 
corresponding to the left-hand side of Eq. (44) can be decomposed into blocks 
of 2 X 2-matrices, which we denote as D(g) for the diagonal part acting on the 
subspacc with wave vector q, and — /pA for the parts which couple subspace 
q with subspace q + 2pQ: 



D(?) = 



2C - 5 -2Ccosq-3fo 



—U) 



\ 



-UJ 



A 



2C -S -2Ccosq- foj 
The eigenvalues of the diagonal part D(g) are given by 

E^{q;u;) ^ 2C - S - 2Ccosq - 2/o T ^u;^ + fl 




(45) 



which determines the downward shift of the unperturbed eigenvalues EQ±{q; uj) 
(cf. Eq. (41)) caused by the static term /o in the expansion (42) (which has 
the same order of magnitude as fi). The perturbation due to \fp\, \p\ > 1, 
will then act mostly by raising the degeneracies of these eigenvalues. The first 
kind of degeneracy. 



E^{q;u) — E+{27r — q;u}) and E_{q;u) — E_{2'k — q;u}), 



(46) 



is raised when q + 2pQ — —q (mod 2%), or 

q — —pQ (mod tt) , 



(47) 



with p integer. Then, a gap of width proportional to appears between the 

two curves E^{—pQ;uj) and E^{pQ;uj), and an identical gap opens between 
the curves E-{—pQ;uj) and E-{pQ;uj) (Fig. 8 (b)). However, the existence of 
these gaps does not change the number of intersections E±{q;u!) = which 
exist in the linear limit, and thus does not induce any instabilities for the SWs. 
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(a) E ▲ 




Fig. 8. (a) Schcmatical picture of the unperturbed band structure E{uj) = EQ±{q] lo) 
(Eq. 41); (b) Gap opening for the first kind of degeneracy (46) (no instability); (c) 
Gap opening for the second kind of degeneracy (48). An instability is generated if 
the gap overlaps the axis E = {). 



For the second kind of degeneracy, at transverse intersection points, a gap 
with width proportional to |/p| opens if 

E+{q- ujp) = E_{q + 2pQ; cUp) = Ep{LUp), (48) 



and a straightforward calculation yields the perturbed eigenvalues at Up to 
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order 



E'^^M = E,{u,) ± ^if=%. (49) 

^p + Jo 



If £■ = belongs to one of these gaps, the SW becomes unstable (see Fig. 8 
(c)). These intersections occur when 



1 

C 



T^V^P + /o =cos(g + 2pg) - cosg (50) 
Ep = 2C - 5 -2fQ-C cos (g + 2pQ) - C cos q. (51) 



It comes out readily from Eqs. (50) and (51) that for a given p, these inter- 
section points are located on an ellipse in the (u;, -E)-plane given by 

^ {2C-5- 2/o - E,f _ 



4C2sin2pQ-/2 4C2cos2pg - /2cot2j9g 



which exists for |/o| < 2C \ sinpQI, with its center at u; = and E = 2C — 
6 — 2/o. Thus, this ellipse represents the locus of the middle of the gaps which 
are opened by the perturbation of order fp, and an instability of the SW is 
generated when this ellipse intersects the axis E = 0, that is when 

{2C-S- 2foY < AC^cos^pQ - f^cot^pQ. (53) 



The existence of these ellipses of gaps is confirmed by numerical calculations 
of the set of bands E{uj) (see Fig. 9). In the incommensurate case, \p\ may run 
from 1 to +00, and thus there are infinitely many ellipses which cover densely 
the rhombic area where the set of curves E-{q; uj) intersects the set of curves 
E_^_{q'] u) for arbitrary q and q'. When E = crosses this area, there are surely 
infinitely many ellipses intersecting this axis, but the instabilities which are 
generated at these intersections become very weak and numerically invisible 
when p becomes large. 

As p = 1 corresponds to the ellipse with the largest gaps, it yields the strongest 
instabilities if these gaps are opening around E = 0. For p = 1 and in the 
limit of zero amplitude SWs (/o = 0,6 = 6o{Q)), the particular ellipse (52) for 



Here, the possible interference with resonances of the first kind (47) has been 
neglected. The case of simultaneous resonances, which occurs at special points, is 
discussed in the Appendix. 
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-;3 -4 -a -1 1 ^ 3 4 b -5-4-3-2-10 1 2 3 4 

CO 03 

Fig. 9. Band spectrum of (39) for a SW with (a) Q/27r = 17/89 ~ ggII {Q < vr/2), 
6' = 1.23, and (b) Q/27r = 34/89 ~ ctg (Q > vr/2), 6' = 3.4 (C = 1). The visible 
gap openings occur along the ehipse (52) with p = 1. 



the intersections {uji, Ei) has the equation: 

ujI {2C cos Q-E^f 

4C2 sin^ Q ^ 4C2 cos2 Q 



1. (54) 



Thus, it is tangent to the axis £" = at u; = 0, and it hes above E = 
for wave vectors < Q < 7r/2 and below for 7r/2 < Q < vr. When the SW 
amplitude is small but nonzero, we obtain using Eqs. (26) and (43) that, 
to lowest order in e^, the 'midgap' ellipse (52) with p = 1 cuts the axis a; = 
at E = — e^/4. Then, using Eq. (49) to calculate the perturbed eigenvalues 
for small but nonzero |ci;i| and |/o| ^ \uji\ we obtain, again by using (43), 
that E[_ ^ — e^/2 and E[_^_ to order when |a;i| is small but nonzero. 
Thus, we should expect that when < Q < 7r/2 (so that the unperturbed 
ellipse is above the axis E = 0), gaps opening in the lower part of the ellipse 
will surround the axis E = for some interval of Ic^l close to uj = O.p^ 
As a consequence, SWs with < Q < n/2 become unstable through first- 
order oscillatory instabilities for arbitrarily small amplitude (Fig. 9 (a)). 
By contrast, as for 7r/2 < Q < tt the unperturbed ellipse lies below the axis 
E = 0, we expect all first-order gaps outside the immediate neighbourhood of 
cij = to open strictly below E = 0. Moreover, as for analytic SWs there are two 
explicit eigensolutions with E' = at = (the phase mode and the sliding 
mode, see Appendix A. 2), there is no gap around E = (and consequently 
no instability) for ti; = 0. We also confirmed by numerical calculations (Figs. 
9 (b) and 10 (a)) that no first-order gaps open around E = 0, so that for 
SWs with 7[ /2 < Q < n no first-order instabilities develop. However, it is 

A similar situation with first order instabihty occurs also for all commensurate 
SWs with Q ^ TT, see Appendix. 
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important to remark that the above estimate based on Eq. (49) is not vahd in 
the immediate neighbourhood of = 0. Actually, in the zero- amplitude limit 
= 0, the point Ei = 0,uJi = where the unperturbed first-order ellipse 
(54) is tangent to the axis £■ = is such that Eq^{Q) = Eq_{0) = 0, which 
imphes by Eqs. (41) that q — ±Q so that also the first resonance condition 
(47) is fulfilled for \p\ — 1. Thus, close to the origin E = 0,uj = the two 
different kinds of gap openings strongly interfere with each other at the same 
order (p = 1), and thus must be considered simultaneously as described in the 
Appendix. 

0.01 
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Fig. 10. Part of the band spectrum of (39) for a SW with Q/2n = 27/89 {Q > 7r/2) 
and 6' = 2.62. The thick dashed line in the left figure is the second-order ellipse 
{p = 2 in Eq. (52)). Right figure shows a magnification of the part around E = 
where a second-order gap opens. 



Thus, when t:/2 < Q < tt, the strongest instability appears for the smallest 
p > 1 such that the condition (53), which in the zero-amplitude limit reduces 
to I cospQl > I cosQl, is fulfilled. An example with p = 2 is illustrated in Fig. 
10. Then, there always exists a p such that, for arbitrarily small amplitude, 
the SW is unstable to order p. Anyway, since the width of the gap openings, 
and therefore also the maximum instability growth rates, are proportional to 
the instabilities become very weak for Q close to tt as the smallest p 
yielding an instability becomes large, and they could be practically invisible 
in numerical simulations for finite-size systems. We might therefore describe 
these analytic SWs as 'quasi-stable'. 



5 Numerical analysis of SWs in KG chains 



Let us now return to the original KG model (5), for which we will investigate 
numerically the validity of our previous results when the conditions of small 
amplitude and small coupling are not necessarily fulfilled. In general, when 
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the SW amplitude is non-negligible also harmonics a^^ with \p\ 7^ 1 in the 
Fourier expansion (10) will affect the dynamics. Moreover, as already men- 
tioned in Sec. 2.2 in connection with Eqs. (15)-(18), when the coupling Ck 
is increased also the small-amplitude dynamics will in general be affected as 
the long-range nonlinear interaction terms in the extended DNLS equation 
become non-negligible, and in particular resonances may occur also for small- 
amplitude SWs if higher harmonics of the solution enter the linear phonon 
band. We will consider the particular cases of the soft Morse potential (3) 
and the hard quartic potential (4). As described in Sec. 3.2, the SWs are ob- 
tained by numerical continuation of the multibreathers with coding sequences 
defined by Eqs. (28)-(29) (modified according to footnote 5 in the case of a 
hard potential), using a modified Newton scheme (described in detail in Ref. 
[11]) to find numerically exact time-reversible solutions {?i„(t)} with a given 
time period Ti, = 2n/u!b. The linear stability of these solutions is then deter- 
mined by numerical integration of the Hill equations (6) over one period Th, 
defining a linear symplectic 2A'^-dimensional map given by the Floquet matrix 



As To is symplectic, for each eigenvalue A there are also eigenvalues at A* 
and 1/A. Then, the solution {un(t)} is linearly stable if and only if all the 
eigenvalues of To({'U„}) lie on the unit circle. The Krein signature associated 
to a complex conjugated pair of eigenvalues A = re='='^ is defined as [32,4] 



where {e„, e„} is the eigenvector of Tq with eigenvalue r e+'^. As before, insta- 
bilities can only occur through collisions of eigenvalues with opposite Krein 
signature. 

At the anticontinuous limit Ck — 0, each site with code |cr„| = 1 yields a 
degenerate pair of eigenvalues at -|-1, while each site with code cr„ = yields 
a complex conjugated pair of eigenvalues at c'^^'^^/^t^ When Ck is increased 
from zero, the degeneracies are raised and the spreading of the eigenvalues is 
just analogous to the DNLS case discussed in the previous section. For type 
H SWs all eigenvalues will initially remain on the unit circle (constituting 
bands for commensurate SWs and Cantor-spectra for incommensurate SWs 
as for the DNLS case), preserving the linear stability until the first oscillatory 
instability appears through a collision on the unit circle between eigenvalues 
of different Krein signature originating from and 6='='^'^/'^'', respectively. 



To(K}): 




(55) 




(56) 
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For type E SWs, in addition pairs of eigenvalues originating from +1 will 
immediately go out on the real axis when Ck becomes nonzero. Thus, in 
general, the instability scenario we observe is as expected identical to that of 
the DNLS model, as long as the amplitude and coupling remain small, and 
no higher-order resonances occur. The scenario outside this regime is however 
quite different for the hard quartic and the soft Morse potentials, and we will 
therefore discuss them separately below. 



5.1 The quartic potential 



Let us first discuss SWs with wave vector Q in the hard quartic potential 
(4), having frequencies LUb > uJi{Q) > 1 (c.f. Eq. (8) and Fig. 5). This case is 
particularly simple since the on-site potential is spatially symmetric, V{u) — 
V{—u), and consequently the time-periodic solutions with fundamental period 
Tf, = 27i/ujh fulfill Unit + Tb/2) = —Unit) (i.c., the time-reversible solutions are 
also time antisymmetric with respect to t = Tf,/4), so that we can restrict our 
numerical Newton algorithm to search for solutions fulfilling this condition. 
Note that this condition is not fulfilled for any even higher harmonic of the 
solution (i.e., with period Tf,/p, p even), so that there can be no resonances 
with even harmonics disturbing the unique continuation of the SW. Actually, 
we already remarked in Sec. 2.2 that there is no coupling between even and 
odd harmonics of the solutions for the quartic potential. Moreover, as V"{un) 
in the Hill equations (6) in this case has the fundamental period Tb/2 rather 
than Tb, it is convenient to redefine the Floquet matrix from (55) as 



{-en{n/2)} 



\{-e^{T,l2)] 



T' 

J- n 



'K(O)}' 

V{^n(0)}y 



(57) 



Then, the eigenvalues of Tq corresponding to sites with code cr„ = will at 
the anticontinuous limit be located at e^™/'^'', i.e., at half the angle of the 
corresponding eigenvalues of Tq. 

Considering now small- amplitude solutions with increasing coupling Ck-, we 
find for different wave vectors Q three qualitatively different scenarios, de- 
pending on whether the second and/or third harmonic of the SW frequency 
uj}j enters the linear phonon band. In general, the condition that the pth har- 
monic of a SW with wave vector Q has entered the phonon band can be written 
as 

^Ji + ACkSiu^^ <Ub< \]\^^Ck. (58) 
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Fig. 11. Eigenvalues of the Floquet matrix Tq for a type H SW with wave vector 
Q = 7r/4 in a finite-size KG lattice {N = 120) with the quartic potential (4), 
ujb = 1-8, and increasing Ck- (a) Ck = 0.199 (stable wave); (b) Ck = 0.699 
(oscillatory instabilities); (c) Ck = 2.773 (additional non-oscillatory instabilities 
due to collisions at —1 resulting from second-harmonic resonances, see inset); and 
(d) Ck = 3.399 (return to stability close to the linear limit for finite-size system). 



This condition is fulfilled only when 

sm — < - and Ck > 
2 p 



P 



4(1 



p2 sin^ ^) 



(59) 



Thus, for SWs with wave vector Q > n/3 there are no higher order resonances. 
In this case, the SWs can for fixed frequency ujb be continued versus Ck until 
they reach their linear limit at Ck — . \ q (where Ub — toi{Q) according to 

4sin -J 

(8)), and the instability scenario during this continuation is completely anal- 
ogous to that observed in the DNLS approximation discussed in the previous 
section, also when ujb is not close to 1. On the other hand, when Q < n/3 
the second harmonic enters the phonon band before the SW reaches its hnear 
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limit if ujh is large enough. As discussed above, this resonance does not affect 
the unique continuation of the SW which still can be continued to its linear 
limit, but causes additional non-oscillatory instabilities which are signalled by 
collisions of eigenvalues of Tq at —1. An example of the instability scenario 
observed in this regime is illustrated in Fig. 11. 

The third regime occurs when sin y < |, or Q ^ 0.2167r. Then, if ujh is large 
enough so that the conditions (58)-(59) are fulfilled for p = 3, the SW with 
frequency Uh will generally bifurcate with solutions with frequency 3co'b be- 
fore reaching the linear limit. These bifurcations are signalled by collisions at 
+1 in the Floquet eigenvalue spectrum of Tq. We will here not attempt to 
give any rigorous description of the possible nature of these (sometimes rather 
complicated) bifurcations, but rather make some qualitative general observa- 
tions (more details are given in [35]). A typical example of a scenario resulting 
from such bifurcations is illustrated in Fig. 12. In this case, the continuation 
from the anticontinuous limit of the SW with wave vector Q = 47r/21 and 
Ub = 4.5 yields the first collision at +1 resulting from a third-harmonic res- 
onance at Ck ~ 44.18 (see Fig. 12 (c) and the magnification in (d)). This 
is slightly before the third harmonic enters the linear phonon band (Eq. (58) 
yields Ck = 45.3125), since the wave amplitude is non-negligible and the res- 
onating wave is a nonlinear rather than a linear SW. As illustrated in Fig. 
12 (a) and (d), it is possible to continue smoothly the solution through this 
bifurcation, but the nature of the solution changes drastically. In fact, this 
scenario typically occurs in two steps: at the first bifurcation point the solu- 
tion changes into an intermediate (unstable) solution, and later, at a second 
bifurcation point {Ck ~ 44.8 in Fig. 12) into the final solution. The latter is 
the SW whose wave vector Q' (close to vr) coincides with the wave vector of 
the resonating eigenmodeP^, and it corresponds to an anticontinuous coding 
sequence generated by Eqs. (28)-(29) but with the coding sequence dn taking 
the values ±3 instead of ±1. Then, the smooth continuation versus Ck for 
fixed Ub follows this new solution without further bifurcations, until it reaches 
its linear limit at Ck = — V?^- 

However, it is important to note that there generally exists also a solution 
with the original wave vector Q close to the linear limit (in smooth continu- 
ation of the linear SW with wave vector Q), but due to the third-harmonic 
resonances there should in this case be no continuous path in parameter space 
connecting this linear SW with the SW with wave vector Q generated from 
the anticontinuous limit. However, if the numerical continuation is performed 
with less care, the resonances might be missed and an apparently (but falsely!) 
smooth continuation of the SW with wave vector Q from the linear to the an- 
ticontinuous limit can be observed as illustrated in Fig. 12 (b)-(c). The more 

^■^ A necessary resonance condition for the quartic model is Q' = mQ mod 2tt, m 
odd. The strength of the resonance generally decreases with |m|. 
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Fig. 12. Bifurcation scenario resulting from third-harmonic resonances in the quartic 
model. The initial type H SW has wave vector Q = 47r/21 and frequency Ub = 4.5. 
(a) Smooth continuation for increasing Ck of a SW with wave vector Q generated 
from the anticontinuous limit into a type H SW with frequency Sw^ and wave 
vector Q' = 207r/21. (b) Intentionally 'careless' continuation from the linear limit 
for decreasing Ck of a SW with wave vector Q, where the step size in Ck has been 
chosen large enough to jump over all resonance regimes (note the small visible gap at 
Ck ~ 44.2 corresponding to the resonance in (a)), (c) Arguments of all eigenvalues 
re'^ of the Floquet matrix Tq versus coupling for one unit cell {N = 21) of the SW 
with wave vector Q in (b) . (d) Magnification of the first resonance in (c) . The solid 
line represents the smooth continuation from below as in (a), the dotted line the 
continuation from above as in (b) which stops at Ck ~ 44.22. (e), (f) Magnification 
of the second and third resonances, respectively (see text for explanation). 
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careful numerics yields the behaviour close to the three resonances illustrated 
in Fig. 12 (d)-(f). Thus, we find that there exists a SW with the original wave 
vector Q also immediately after the first resonance, but attempts to continue 
it towards smaller Ck (dotted line in Fig. 12 (d)) fail for Ck ~ 44.22. Sim- 
ilarly, the attempts to continue it towards larger Ck fail when approaching 
the second resonance at Ck ~ 46.51932 (left part of Fig. 12 (e))P^, so that 
its existence region apparently is limited by these two resonances. A similar 
behaviour is observed also between the second and the third resonances. Ap- 
proaching the second resonance from above the continuation of the SW with 
wave vector Q stops at Ck ~ 46.51935 (right part of Fig. 12 (e)), while ap- 
proaching the third resonance from below it transforms, via an intermediate 
solution appearing through a collision at -|-1 at Ck ~ 51.826293472 (solid line 
in Fig. 12 (f)), into a SW with wave vector Q'" = 167r/21 similarly as in Fig. 
12 (a), (d). Finally, the smooth continuation of the linear SW with wave vec- 
tor Q towards smaller Ck stops when approaching the third resonance from 
above (open circles in Fig. 12 (f)) at Ck ^ 51.826293474688. (As the third 
resonance is extremely weak, its correct numerical resolution requires working 
with quadruple precision and using a step size in Ck smaller than 10~^^.) 



Let us finally note that in all cases, we find that for type H SWs in finite 
systems with periodic boundary conditions, stability will always be recovered 
in an interval close to the linear limit (c.f. Fig. 11 (d)), which shrinks to zero 
as the system size increases. The return to stability occurs through reentrant 
instabilities, completely analogous to the scenario for the DNLS approximation 
mentioned in Sec. 4.1. The shrinking of the stability regime with system size 
is illustrated in Fig. 13. 

5.2 The Morse potential 

For the soft Morse potential, the nonlinear SWs with wave vector Q have 
frequencies Ub < u)i{Q), and the continuation from the anticontinuous to the 
linear limit must be performed in two steps as illustrated in Fig. 5: first by 
increasing Ck for fixed < 1, and then by increasing ujf, for fixed Ck- Even 
though the on-site potential (3) now is asymmetric, we can use the symmetry 

However, a smooth (but non-monotonous) continuation into a SW with wave 
vector Q" = 187r/21, analogous to that in Fig. 12 (a), can be obtained in the 
following way [35]: at Ck ~ 46.51932, there is a bifurcation with a second solution, 
which can be followed smoothly by decreasing Ck until (at Ck ~ 46.51914) it 
bifurcates with a third solution, which, again by increasing Ck, can be followed as 
it smoothly changes into the SW with wave vector Q" . This 'Z-type' behaviour, with 
the simultaneous existence of three different solutions for the same Wf, and Ck in 
some regime, typically occurs in systems with effective ultra-long-range interactions 
yielding length-scale competition [36]. 
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Fig. 13. Lower limit for the interval of stability close to the linear limit as a function 
of system size, for a type H SW with Q = 47r/7 and ojf, = 3.0 in a quartic KG 
lattice. The dashed line is the linear limit according to Eq. (8). 

properties of the SW coding sequences discussed in Sec. 3.2 to restrict the 
space of solutions for the numerical Newton scheme. To find the time-reversible 
SWs of type H with coding sequences antisymmetric around a lattice site 
we impose the conditions M„(Tfc/2) = m_„(0), Un{Th/2) = Un{0) = in the 
Newton algorithm, while for SWs antisymmetric around a bond center the 
first condition is replaced by Un{Tb/2) = u-n+i{0). For type E SWs, the 
symmetry condition Un{t) — U-n+i{t) can always be employed. 



In the case of a soft potential, the pth harmonic of a SW with wave vector Q 
enters the linear phonon band if (c.f. (58)) 



- < a;^ < mm 
P 



(60) 



Thus, in contrast to the hard case, we expect to find higher order resonances 
for all Q if, for fixed Ck, the frequency uj^ is small enough, i.e., if the SW 
is sufficiently nonlinear. In particular, when the conditions (59) are fulfilled, 
the resonances occur immediately from the linear limit. As the potential is 
non-symmetric already the second-harmonic resonances yield generally bifur- 
cations of the SWs, and thus we find qualitatively different behaviour for 
small-amplitude SWs depending on whether the condition (59) with p = 2, 

TT 3 

Q<- and Ck>-j- ^ . , , (61) 

3 All- A sm^ f j 
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Fig. 14. Phase diagram for SWs with Q = 47r/5 in the Morse KG model. The lower 
solid line is the linear dispersion curve (8), the filled circles indicate the line (which 
in fact are two closely separated lines, see text) of inversion of stability between the 
SWs of types E and H (characterized by eigenvalue collisions of To at +1), and the 
open squares indicate the line of amplitude divergence. The upper solid line is the 
right boundary of the second-harmonic band. Inset shows the difference between 
the values of the coupling where the stability changes for the H and E waves for 
fixed frequency. 

is fulfilled or not. In order to avoid complications resulting from higher-order 
resonances {p > 2), which occur only for highly nonlinear SWs, we will mainly 
concentrate the discussion here to the regime uJb > 1/2 (more details concern- 
ing the behaviour for smaller Ub are given in [35]). 

The typical behaviour when continuing SWs with Q > 7r/3 is summarized 
in Figs. 14-15, for a rather simple case with Q = An/ 5. Let us first concen- 
trate on the part closest to the linear limit. Then, we find that increasing 
Ck for a fixed uJb close to 1 (which means that the wave amplitude also in- 
creases), the type H SWs will suffer additional real instabilities, not present 
neither in the DNLS approximation nor in the quartic model. These instabili- 
ties occur through collisions at +1 at some critical value of the coupling. Very 
close to this value, also the type E SW with the same Q and Ub changes its 

In fact, the ansatz (28)-(29) for the SW coding sequences is valid only for 
LJb > 1/2; a smooth continuation of SWs generated from (28)-(29) for tOb > 1/2 
towards smaller ivb shows [35] that they appear with a different anticontinuous cod- 
ing sequence when ujb < 1/2, with the O's in the original code replaced by ib2's. 
Solutions generated from the anticontinuous limit at Wfc < 1/2 using (28)-(29) are 
lost through bifurcations and not continuable to the linear limit [35]. 
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Fig. 15. Continuation from the linear limit at fixed coupling Ck = 0.8 for decreasing 
frequency for (a),(c) the type H and (b),(d) the type E SW with Q = 47r/5 in the 
Morse KG model (corresponding to a horizontal line in the phase diagram in Fig. 
14). The wave amplitude diverges in both cases at ujh ~ 1.701. (c) and (d) show the 
magnitude of the eigenvalues re*^ of Tq corresponding to one unit cell (A'^ = 5) of 
the H and E waves, respectively. Insets are magnifications of the region between 
the linear limit and the point of inversion of stability, which occurs at ~ 1.9155. 
The other instabilities are oscillatory instabilities resulting from Krein collisions. 

stability properties, in that the real eigenvalues which are normally present 
collide at +1 and go out on the unit circle, leaving only Krein instabilities and 
second-harmonic collisions at -1 as possible remaining instability mechanisms. 
This inversion of stability' is analogous to the scenario occurring between the 
single-site (code ...0001000...) and two-site (code ...0011000...) breathers in the 
Morse and other soft asymmetric potentials [37,38], and is generally associ- 
ated with a highly increased mobility [38,20,21]. In the case shown in Figs. 
14-15, the change of stability for the ii^-wave with a given frequency occurs 
for slightly smaller values of Ck than that of the corresponding H-waye (see 
inset in Fig. 14), so that in the absence of other instability mechanisms (which 
is the case for small systems close to the linear limit, c.f. Fig. 15), there is a 
small regime in the {ub, Cii-)-plane of simultaneous stability for both wavesf^. 

By contrast, for the single breather there is a small regime of simultaneous insta- 
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When increasing Ub and Ck the regime of simultaneous stabihty decreases, 
and moreover at a certain point {uf, ~ 2.1 and Ck ~ 1-0 for the case in 
Fig. 14) the two hues defining the inversion of stabihty converge and intersect 
the hnear dispersion curve (8). Thus, for values of Ck above this point the 
stability scenario close to the linear limit is opposite to what we have found 
previously in the DNLS-approximation: the H-wave becomes immediately un- 
stable through real eigenvalues, while the ii^-wave remains stable in an interval 
close to the linear limit for finite-size systems. 

When further increasing the coupling, we find typically that the region of ex- 
istence for the SWs is limited by a curve in the {uh, Ci^)-plane where the wave 
amplitude diverges as illustrated in Fig. 15 (a), (b) (as the Morse potential (3) 
becomes flat for large u). However, for smaller Ub the scenario becomes more 
complicated with additional bifurcations, and we will not attempt to give an 
exhaustive analysis of this regime here (see [35] for more details). Let us just 
mention that for the case illustrated in Fig. 14, when Ub < 0.72 {u^ < 1.1) for 
the if- wave (ii^-wave), the continuation from the anticontinuous limit at fixed 
Ub is interruptedp^ by a bifurcation before reaching the point of amplitude 
divergence. Let us also mention that for SWs with more complicated unit cells 
than that of Figs. 14-15 (which is somewhat special since the coding sequence 
for one unit cell of the if-wave [1-101-1] has only one per unit cell, while 
that of the E-wave [-11-11-1] has no and only one pair of consecutive -1), 
there are sometimes more than one 'stability-inverting' bifurcation, and the 
phase diagram of 'stable' (neglecting other instability mechanisms) and un- 
stable regimes in the {ub, Ci^)-plane may contain more than two parts [35] (as 
was also found for the single breather [38]). 

Let us now turn to the case with Q < 7r/3, in which case we find second- 
harmonic resonances also in the linear regime for large enough coupling. The 
scenario in the small-amplitude regime is then similar to the scenario for the 
third-harmonic resonances in the quartic model illustrated in Fig. 12. An ex- 
ample for a type H SW with Q = 27r/ll is shown in Fig. 16. Continuing from 
the anticontinuous limit for increasing coupling then yields (Fig. 16 (b)), via 
an eigenvalue collision at +1 for a value of Ck on the line of open diamonds in 
(a), a transition to an intermediate state, which at a second bifurcation point 
becomes a SW with frequency 2uJb and wave vector Q', where Q' is the wave 
vector closest to vr giving a SW whose symmetry is compatible with that of 
the original wave (here, the new wave is the type E SW with Q' = 107r/ll)P^. 

bility [38]. 

However, it can be restored [35] by following a curvilinear, 'Z-shaped' path simi- 
larly as described for the quartic model in footnote 13. This behaviour is likely to be 
related to a strong effective long-range interaction in the extended DNLS equation 
(12); cf. [36]. 

^"^ Similarly to the quartic model (footnote 12), a necessary resonance condition is 
Q' = mQ mod 27r, but with no restriction on the integer m. 
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Fig. 16. Second-harmonic resonances for a type H SW with Q = 27r/ll in the 

Morse KG model, (a) shows the bifurcation hnes for the first two resonances; dashed 
hnes are the hnear second-harmonic dispersion curves for Q = lOvr/ll (lower) and 
Q = Stt/II (upper); fat solid line is the fundamental linear dispersion curve for 
Q = 27r/ll. (b) illustrates the smooth continuation versus coupling from the anti- 
continuous limit at fixed frequency cj;, = 0.9 through the first resonance (open dia- 
monds in (a)) into a type E SW with frequency 2iOh and wave vector Q' = IOtt/II. 
(c) shows smooth continuation versus frequency from the linear limit at constant 
coupling Ck = 1.4 (i.e., between the first and second resonances) through the sec- 
ond resonance (filled circles in (a)) into a type E SW with frequency 2ujb and wave 
vector Q" = Stt/II. (d) shows intentionally careless continuation which jumps over 
the resonance in (c); if a careful continuation versus larger uob is attempted, it stops 
at LOb 1.18235 (filled circle in (a)) similarly as in Fig. 12 (d) and (e). 

This new SW can then be continued for larger coupling until its amplitude 
diverges at some finite value of Ck (similarly as in Fig. 15). 

The same scenario is observed by continuation from the linear limit for fixed 
coupling towards lower frequency, as long as Ck is below the value where 
the first resonance curve (which in the linear limit coincides with the linear 
second-harmonic dispersion curve for Q = Q', see Fig. 16 (a)) intersects the 
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linear dispersion curve. For Ck larger than this value the SWs are also continu- 
able from the linear limit towards lower frequency, but attempts to continue 
these waves to the anticontinuous limit by subsequently decreasing Ck fail 
since, similarly as for the quartic model (Fig. 12 (d), (e)), the continuation is 
interrupted when reaching the first resonance curve. Performing instead the 
continuation for fixed coupling, the scenario is completely analogous to that of 
smaller Ck'- as illustrated in Fig. 16 (c) for Ck = 1-4, the wave bifurcates at a 
point on the line of filled circles in (a) (0;^ ~ 1.1816) into an intermediate state, 
which in a second bifurcation (at a;^ ~ 1.1789) turns into the second-harmonic 
type E SW with wave vector Q" = 87r/ll, which then can be continued until 
its amplitude becomes infinite. Similarly, if the continuation from the linear 
limit is performed for Ck above the intersection point with this second res- 
onance curve, the continuation towards smaller coupling will be interrupted 
by the second resonance curve, while the continuation towards larger coupling 
or smaller frequency meets a third resonance curve where the wave changes 
smoothly into a SW with wave vector Q'" — 67r/ll, etc. However, as illus- 
trated by Fig. 16 (d), these resonances may be jumped and an apparently 
smooth continuation of the SW with the original wave vector Q observed if 
the numerical continuation is performed taking large steps. 

Finally, we remark that also in the case of second-harmonic resonances close 
to the linear limit, the inversion of stability between the E and H waves 
can be observed. For example, in the case considered in Fig. 16 (c), the first 
eigenvalue collision at -|-1 in the continuation from the linear limit occurs 
at uji, 1.1822 (i.e., very close to the coUision corresponding to the second- 
harmonic resonance), and corresponds to inversion of stability. Let us also 
remark that, as mentioned above, the bifurcation scenario for smaller ujh be- 
comes more complicated, and for example the monotonous continuation for 
fixed frequency versus coupling from the anticontinuous limit of the SW in 
Fig. 16 is interrupted before reaching the first second-harmonic resonance for 
< 0.7. 



6 Dynamics 

Let us now discuss the dynamics resulting from the oscillatory instabilities 
of the SWs over different time ranges. We will here emphasize generic rather 
than model-dependent properties, and therefore we consider KG systems in 
the regime of small amplitude and small coupling, where the DNLS approxi- 
mation is expected to be good and the model-dependent phenomena discussed 
in the previous section (i.e., higher-harmonic resonances and inversion of sta- 
bility) do not play any role. As was illustrated already by Fig. 4 in [25], 
and is further exemplified by the figures in this Section, the dynamics in the 
regime of oscillatory instability can then typically be divided into roughly 
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Fig. 17. Spontaneous motion (resulting only from numerical truncation errors) of a 
single discommensuration in a Morse potential KG chain (5) with Ck = 0.03. The 
initial frequency is uj^ = 1.036. 

three regimes. Initially, the small perturbation grows exponentially leading to 
oscillatory dynamics with small but increasing amplitude. As the oscillation 
amplitude reaches some threshold, an intermediate regime appears which typ- 
ically is characterized by an inhomogeneous translational motion, where the 
wave remains locally coherent but with its different parts moving with respect 
to each other. Gradually, the wave loses coherence and a final, apparently 
chaotic state appears, the statistical properties of which depending crucially 
on the SW wave vector as discussed below. 

To understand the appearance of the intermediate regime of translational 
motion, it is useful to first consider the case Q —>■ tc with one single discom- 
mensuration in the Q = n phonon (for a soft potential). As was found in 
[18,22] for the DNLS model, the oscillatory instability leads in this case to 
the motion of the discommensuration (for the DNLS model interpreted as 
a discrete 'grey' soliton) with a rather well-defined velocity for long times, 
although it continuously emits radiation and decays. The same type of be- 
haviour occurs also in KG models for small Ck as illustrated in Fig. 17, 
and is basically a consequence of the fact that the unstable eigenvector is 
spatially antisymmetric with respect to the background Q = n wave around 
the discommensuration site. Thus, if we regard a SW with Q close to vr as 



In fact, it is possible [39,35] to calculate numerically exact solutions corresponding 
to moving discommensurations with given velocities, although they typically acquire 
an oscillatory extended tail and could therefore be described as 'dark nanopterons' 
(for the DNLS model, moving dark solitary waves close to the continuum limit have 
been numerically calculated in [40]). 
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Fig. 18. Time evolution over short and intermediate time scales of (parts of) slightly 

perturbed SWs with Q/2tt = 34/89 ~ ac in a Morse potential KG chain (5) with 
(a),(c) Ck = 0.03 and = 1.034 (nonanalytic case, type H), and (b),(d) Ck = 0.05 
and LOb = 1.075 (analytic case), (c) and (d) illustrate the time evolution and inter- 
action of 'discommensurations', defined as sites where the energy is (c) Hn < 0.01 
respectively (d) Hn < 0.00232 (corresponding to sites with code for the initial 
SWs). Note the very slow evolution of the instability in the analytic case, which 
can be considered as 'quasistable' as discussed in Sec. 4.2: the strongest instability 
appears only to order p = 4 according to Eq. (53) . 



an array of interacting discommensurations (where each discommensuration 
corresponds to a site with code 0), the intermediate regime of mainly transla- 
tional motion can be anticipated. As illustrated in Fig. 18, this description is 
at least qualitatively relevant also for analytic waves (where generally discom- 
mensurations can no more be regarded as localized objects as their range of 
interaction diverges [19]) with wave vectors not necessarily close to tt (so that 
the discommensuration sites may be close). Initially, the collective movement 
of the discommensurations is well-ordered and they remain clearly separated 
from each other, but as time evolves their motion becomes more and more 
erratic and they start to collide and possibly merge with each other. Thus, 
the SW eventually loses its coherence. 
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Fig. 19. Time evolution over short and intermediate time scales of (parts of) a 
slightly perturbed SW with Q = 127r/55 (< vr/2) and = 1.0015 (corresponding 
to an analytic wave in the incommensurate limit) in a Morse potential KG chain (5) 
with Ck = 0.05. (b) illustrates the time evolution and interaction of high-energy 
sites ('breathers'), defined as sites where the energy is Hn > 0.012 (corresponding 
to sites with codes ±1 for the initial SW). 

When (for a soft potential) Q < 7i/2, the coding sequences contain consecu- 
tive zeros, and therefore the discommensurations in the initial SWs are not 
isolated. Thus, for small Q it is more natural to regard the SWs as composed 
of individual breathers, where each breather corresponds to a site with code 
±1, and the distance between the breathers increases as Q decreases. Then, as 
illustrated in Fig. 19, the intermediate regime can be interpreted as consisting 
of moving breathers which collide inelastically with each other. 

Let us now discuss the nature of the asymptotic long-time states for different 
Q. For the DNLS equation (19), which is non-generic among nonlinear lattice 
equations in the sense that it has in addition to the Hamiltonian 

n {{^n} , {-i<}) - E -^\^n\' + TjlV^nl' + C\i^n+1 - ^nH (62) 
n=l ^ ^ ^ 

also a second conserved quantity, the excitation norm 

N 

J\f=Y.M', (63) 

n=l 

it has recently been shown [23] that the nature of the asymptotic dynamics for 
typical initial conditions depends critically on the values of these two conserved 
quantities. Redefining the Hamiltonian (62) as H' = 7i+{S—2C)N' (equivalent 
to applying the gauge transformation ipn — > "^n 6"'^''"^*^^* to the conjugated 
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variables {ipn} , {~^4'n})y ^ phase transition was found to occur along the line 



in the thermodynamic limit N oo. On one side of this line (below for a > 
and above for cr < 0), the system thermalizes in a standard Gibbsian sense with 
well-defined temperature and chemical potential, and the equilibrium proba- 
bility distribution function (PDF) for l^ipnl^ decays rapidly for large lipnl"^ so 
that large-amplitude excitations are highly improbable. However, on the other 
side of the line (64) the system exhibits a negative-temperature type behaviour 
where the stationary PDF for {ipnl"^ shows a positive curvature, corresponding 
to the creation of persistent localized high-amplitude standing breathers. The 
phase transition line corresponds to the limit of infinite temperature in the 
Gibbsian regime, and was obtained by assuming independent thermalization 
at each site, i.e., by neglecting the coupling term in the Hamiltonian. Now, we 
note that for a type H SW with Q = 7^/2 the coupling term is exactly zero, 
since according to the code sequence (29) such a wave must have the struc- 
ture 4>2m+i = 4'2m = for all 5. Thus, the Q = tt/2 SW lies exactly 
on the phase transition line (64). Moreover, as illustrated in Fig. 20, this im- 
plies that type H SWs with a given Q always belong to the same well-defined 
phase regardless on the value of 6', and that SWs with Q < 71/2 and those 
with Q > 71 /2 always belong to different phases.\^ By regarding a SW with 
wave vector Q = 7\r' j s' as a multibreather with r' breathers (or, equivalently, 
s' — r' discommensurations) in each (half-)period of s' sites, we obtain close 
to the anticontinuous limit (i.e., asymptotically for large A/"/A^):P^ 



where 9(x) is the usual step function. 

Thus, within the DNLS approximation and assuming a soft potential [a = 
— 1), we should expect for large systems to always observe an equipartitioned 
asymptotic state for all initially slightly perturbed unstable type H SWs with 
Q > 71 /2, while persistent breathers should be created from SWs with Q < 
7r/2 (with the opposite conclusion for hard potentials). We confirmed this 
by numerical simulations, although for Q close to but smaller than 7c/2 the 
system size N needs to be rather large to allow breather formation. Moreover, 
we confirmed numerically that also for the original KG lattices there appears to 
be a phase transition at Q = 7r/2 between a thermalized phase and a 'breather 

Note that this is not in general true for propagating waves, which for a given 
wave vector Q can cross the transition line (64) as the amplitude is varied [23]. 
'^^ This expression is exact for Q = 7r/2 and Q = tt. 
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Fig. 20. Average Hamiltonian 7i' /N versus average norm J\f /N for type H DNLS 
SWs {a = —1, C = 1) having wave vectors (increasing from left to right) Q = vr/S, 
7r/4, 37r/8, 7r/2, Svr/S, 37r/4, 77r/8, and tt. The points corresponding to Q = tt/2 
(filled squares) coincide exactly with the phase transition line (64). 




100000 



Fig. 21. Long-time dynamics of slightly perturbed unstable type H SWs in a 
Morse KG lattice with Ck = 0.035 and = 189. (a) shows breather forma- 
tion for Q = 27r/7 < it/2 and = 0.995, while (b) shows thermalization for 
Q = 6tt/7 > 7r/2 and Ub = 1.035. 



phase' (Fig. 21). However, it should be stressed that the results derived in [23] 
depend crucially on the existence of a second conserved quantity, and it is 
therefore questionable whether such a transition strictly speaking occurs also 
in other models, or if the 'breather phase' observed here for KG models only 
describes the dynamics over a long but finite time range, while asymptotically 
thermalization will occur also there. In any case, over the time scales that 
we have been able to follow numerically we have seen no signs of breather 
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decay, and thus it is clear that the phase transition obtained in the DNLS 
approximation is of large relevance also for the long-time dynamics in the KG 
models. 



7 Conclusions 

In conclusion, we have shown that for a general class of anharmonic lattices, 
different types of SWs with spatial period commensurate or incommensurate 
with the lattice exist, but are generically unstable through oscillatory instabil- 
ities for small amplitudes. Thus, in discrete systems, the nonlinearity- induced 
coupling between counter-propagating waves leads to wave break-down for 
arbitrarily weak anharmonicity even if the individual propagating waves are 
modulationally stable. The physical mechanisms of the oscillatory instabilities 
can be related to a two-component structure of the SWs, which however is of 
different nature close to the anticontinuous and to the linear limits, respec- 
tively. In the first case, the natural interpretation is in terms of the multi- 
breather description (28)-(29) of the SW, with (quasi-)periodically repeated 
sites of large (code ±1) and small (code 0) amplitude, which defines a two- 
sublattice structure. Then, for small inter-site coupling the eigenmodes of the 
linearized equations describing the small- amplitude oscillations around the 
exact SW are mainly localized on one of these two sublattices, and the insta- 
bilities appear when, as the coupling is increased, eigenfrequencies belonging 
to the different sublattices start to coincide. On the other hand, close to the 
linear limit, the oscillatory instabilities arise from overlap between positive- 
and negative-frequency plane-wave solutions to the linearized equations (32) 
in the DNLS-approximation, which are coupled by the SW amplitude tp"^. 

It is instructive to point out some aspects where the stability analysis of SWs 
considerably differs from that of propagating waves. Firstly, for propagating 
waves there is no two-component structure at the anticontinuous limit, as 
their amplitude is homogeneous with all sites having codes |5"„| = 1 [7]. Sec- 
ondly, propagating waves are non-time-reversible with a complex amplitude 
ipn = e^*^" in the DNLS approximation, which modifies the eigenvalue prob- 
lems (32)-(33) [31]. Thirdly, the spatial Fourier expansion (42) of a propagating 
wave in the DNLS approximation trivially contains only one term {p = 1), so 
that instabilities resulting from higher-order resonances (which are the only 
instabilities appearing for analytic incommensurate DNLS SWs close to the 
linear limit) are excluded. 

The existence of a regime of two-channel linear phonon scattering in the DNLS 
equation was also discussed in Ref. [41]; however, in contrast to what was claimed 
in these papers this regime does not exist for a single localized breather (the overlap 
described by Eq. (35) vanishes as Q ^ 0). 



47 



We should also stress that in general, there is no direct relation between the 
frequency of the SW and the frequencies of the oscillatory instabilities, as the 
instabilities are related to resonances between two different classes of solutions 
to the linearized equations (with opposite Krein signature), and not to the SW 
frequency itself. Also, as the instabilities for the infinite chain result from over- 
lap between continuous bands and not isolated eigenvalues, the wavelengths of 
the unstable modes generally belong to a broad spectrum. However, for com- 
mensurate waves close to the linear limit, we proved in the Appendix that the 
first-order (strongest) instabihties correspond to frequencies and wavelengths 
close to, but different from, those of the SW itself {to ^ 0, |g| fa Q, cf. Eq. 
(A. 26)). Moreover, we showed by perturbation theory that the strength of the 
oscillatory instability (i.e., its growth rate) close to the linear limit is pro- 
portional to the corresponding spatial Fourier component \fp\ of its squared 
DNLS amplitude -0^, and thus initially increases with the amphtude. The in- 
stability persists until, in the highly nonlinear regime, the frequencies of the 
two subclasses of linearized modes become out of resonance (through an in- 
verse Krein bifurcation), so that close to the uncoupled limit the instability is 
suppressed (cf. the similar scenario for the discrete dark soliton corresponding 
to Q ^ TT in [18]). 

Obviously, this paper has been mainly concerned with fundamental properties 
of the SWs and their instabilities in generic but simple models, and more work 
needs to be done in order to elucidate their role in real physical systems, e.g., 
whether they have some influence on thermal transport coefficients. However, 
one result obtained here (Sec. 6) that should be emphasized in this context 
is, that the partition of the phase space for the asymptotic dynamics of the 
DNLS model (due to the norm conservation) in two regimes found in Ref. [23], 
corresponding respectively to a Gibbsian thermal equilibrium and formation 
of persisting breathers, is given a clear physical interpretation in terms of the 
SWs: the transition hne (64) is defined by a (type H) SW with wave vector 
Q — 71 / 2, and SWs with given Q always belong to the same regime regard- 
less of their amplitude. More work should be done in order to elucidate the 
relevance of this transition for the KG system, e.g., to see if the partition of 
the phase space is fully broken in the absence of norm conservation, and to 
determine typical hfetimes of 'special' configurations hke breathers and SWs. 
In this context, we also want to draw the attention to the very recent Ref. [42], 
where the processes of breather formation, coalescence and evolution to en- 
ergy equipartition were analyzed for the FPU chain. There, it was found that 
exciting initially a single linear mode, the initial state first relaxes to an anhar- 
monic spatially periodic solution, which we interpret as the FPU counterpart 
of our SWs. Later, this solution breaks up and the periodicity is lost, which 
we believe is related to instabilities similar to those obtained in our work. 
Thus, the SWs and their instabilities apparently provide an important link in 
the understanding of the pathway to energy localization and equipartition in 
anharmonic lattices. 
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Concerning possible experimental observation, the general nature of our results 
suggests that these instabilities should be observable in macroscopic as well 
as microscopic contexts. In particular, as was recently experimentally verified 
[43], nonlinear optical waveguide arrays provide a direct application of the 
DNLS equation, and thus these systems are also good candidates for detection 
of SW instabilities. 
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A Linear stability analysis of small-amplitude standing waves close 
to double resonances 

A.l Points of simultaneous resonances 

Consider the eUipses (52) describing, for each p, the points in the (a;, E)- 
plane where pth order resonances of the second kind (48) occur. On each such 
ellipse, there are special points where also the first resonance condition (47) is 
fulfilled, so that both the above resonances exist simultaneously. These points 
are obtained for special values of q in the parametrized representation of the 
ellipse {ujp{q) , Ep{q)) defined by Eqs. (50) and (51) with parameter q — p'Q 
ip' integer). Then, a gap with a width of the order of \fpi \ opens in the ellipse, 
which makes this perturbed ellipse discontinuous. Thus, even if the ellipse 
(52), describing the center of the pth order gap openings of the second kind, 
intersects the real axis £■ = so that an instabihty should be expected (i.e., 
if condition (53) is fulfilled), this disconnected curve might not intersect the 
real axis if £" = is close to a point of simultaneous resonance, so that there 
is not necessarily an instability. However, since fp decreases exponentially 
as a function of p, only the gaps generated for \p'\ < \p\ makes a significant 
discontinuity in the ellipse of pth order gaps. Thus, as the number of significant 
discontinuities in each ellipse is finite, we expect that only in exceptional cases 
a predicted gap opening around E = Q might be prohibited, and consequently 
the condition (53) for pth order instability should be generic. 

The values of Q for which simultaneous resonances of orders p' and p occur 
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on the axis E = are obtained by vanishing Eq. (51): 

2C-6- 2/0 -Ccos{p + 2p)Q -CcospQ = 0. (A.l) 

In the zero-amphtude hmit fo = 0, 6 = So{Q) (Eq. (24)), this condition reduces 
to 

cos Q = cos pQ cos {p + p')Q, (A. 2) 

which is only fulfilled for a discrete set of values of Q when \p\ > 1. However, 
for the ellipse (54) with p = 1, the condition (A. 2) is fulfilled for all Q when 
p' = —1. Thus, in the zero-amplitude limit there is always a point of double 
first order resonances at = where the corresponding ellipse is tangent to 
the axis E = 0. Consequently, we should treat this case separately as done 
below. 

We will here discuss both commensurate and incommensurate SWs. In the 
commensurate case with Q = Tir'/s', {ipn} s'-periodic sequence, and 2pQ 
modulo 2iT takes only s' different values. Then, the Fourier expansion (42) of 
ip^ will be a finite sum, where \p\ runs only from to s' /2 when s' is even, and 
to (s' — l)/2 when s' is odd. Consequently, there is only a finite number of 
ellipses (52) in the commensurate case. Furthermore, it follows from Eq. (29) 
and the model symmetries that the sequence {ipn} ^ periodic cycle will be 
symmetric either around a lattice site, in which case we choose ipn = i'-ny 
around a bond center, so that we choose ipn — V'-n+i-Q ^^^^ case, we 

have the phase = in Eq. (42), so that fp = f-p are real. In the second 
case, = —Q/2 in (42) so that fp = f*p in general would be complex, but 
by redefining the Fourier transforms (40) and (42) by changing n — n — | we 
again have = and fp = f_p real. 

A. 2 General properties of the hand spectrum at u = 

Let us consider the eigenequation (38) for the special parameter value u = O.P^ 
Then, it splits into two independent eigenequations which can be solved either 
with Wn = or Un = 0: 

C,U^ = {2~6- 3^l)Un - f/„+i - Un-l = EUn (A.3) 
CoWn = {2-6- ijl)Wn - lV„+i - Wn-i = EWn. (A.4) 

^■^ Both symmetries may also exist simultaneously, c.f. the discussion in Sec. 3.2. 

For convenience, in the rest of this Appendix we assume C = 1 and a = —1 
without loss of generality. 
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Eq. (A. 3) can be related to the nonlinear symplectic map (23), for which the 
linearized equation is 

Lxen = (2 - 5 - 3^/'^)en - e„+i - en-i = (A. 5) 

where e„ is real and time constant. As well as for Eq. (30), we can associate 
with this linear equation a real eigenvalue problem defined by £ie„ = Ee^^ 
which is nothing but Eq. (A. 3). It was shown very generally in [44] that a 
necessary and sufficient condition for a trajectory of a symplectic map to 
be uniformly hyperbolic is that = does not belong to the spectrum of the 
operator defined by the second variation of the action, which in our case is L\. 

When ■^^ is space periodic, and thus represented by a periodic cycle in the 
map defined by Eq. (23), the spectrum of Eq. (A. 3) consists of bands separated 
by gaps. Then, if (and only if) the cycle is hyperbolic = belongs to a 
gap, if it is parabolic is at the edge of a gap, while for an elliptic cycle 
belongs to the inner of a band. When {'?/'„} is represented by a cantorus in 
the map (23), the corresponding trajectory is uniformly hyperbolic and then 
does not belong to the spectrum of Eq. (A. 3). On the other hand, when it is 
represented by a KAM torus belongs to the spectrum, and moreover there is 
an explicit eigenvector for this eigenvalue. In that case, ■i/'n is a quasiperiodic 
function of n with an analytic hull function Xs(^) defined by Eq. (27), and 
Eq. (A. 3) exhibits the solution C/„ = dxsiQf^ + 0)/<90 (the shding mode) for 
the eigenvalue E — 

The second eigenequation (A. 4) has always the explicit solution Wn = ipn 
(the phase mode) for the eigenvalue E = 0. When ip"^ is space periodic with 
period s', the eigensolutions of Eq. (A. 4) are Bloch states with a wave vector qgi 
defined modulo 2n/s', and the corresponding eigenvalues form bands separated 
by gaps. Since the solution Wn = ipn has the wave vector qs' = or vr modulo 
l-n j s\ its corresponding eigenvalue = is at the edge of a gap of the 
spectrum of Eq. (A. 4). When {t/'n,} is space quasiperiodic and represented by 
a KAM torus in the map (23), the spectrum of Eq. (A. 4) is a Cantor set with 
infinitely many gaps. Again, = is at an edge of a gap of this Cantor 
spectrum. 



A. 3 First- order perturbative approach for uo and E dose to zero 

We now take advantage of these informations for studying the spectrum of 
Eq. (38) close to the point of double first-order resonances {uJ,E) ^ (0, 0) 
in the limit when the SW amplitude ipn is small and thus close to a cosine 
function. Then, at the lowest order p— 1, it is sufficient to assume that /o 7^ 0, 



51 



/i = /_i 7^ and /p = for |p| > 2. In analogy with Eq. (43) we set 

/o = 2/i + «, (A.6) 



where \a\ <S /o but in general not strictly zero. Then, for describing consis- 
tently at leading order the band structure in the neighbourhood of the origin 
in the (cj, £')-plane, we should essentially consider the reduction of Eq. (44) to 
the subspace involving the four almost degenerate eigenvectors which are the 
two components with wave vector q and the two components with wave vector 
q + 2Q with q ^ —Q. Then, close to the origin uj — 0, E — 0, the bands of 
eigenvalues Ex[q]uj) are well described by diagonalization of the 4 x 4-matrix 



M(g;cu) 



(a{q)-fo -co -3/i 

-u; a{q) + fo -/i 

-3/i a{q + 2Q)-fo -u 



\ 



\ 







-/i 



(A.7) 



a{q + 2Q) + fo J 



where 

a{q) = 2 - 5 - 2/o - 2 cos g. (A.8) 



A. 3.1 Spectrum at u — 

When = 0, we should recover from Eq. (A.7) the properties of the spectrum 
of Eq. (A. 3) associated with the tangent map (A. 5), as well as the phase mode 
of Eq. (A. 4). Then, the 4x4 matrix (A.7) splits into two diagonal blocks of 
2x2 matrices. 



L5(g) 



«(?) - fo -3/i ^ 
-3/i a{q + 2Q)-foj 



(A.9) 



which yields the spectrum of (A. 3) calculated in perturbation for q close to 
—Q as 



Es±{q) = 2 - 5 - 3/o - cos ? - cos {q + 2Q) 
±^{cos q - cos {q + 2Q)Y + Q/f, 



(A.IO) 
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and 



Lp{q) = 



+ fo -fi 
-h a{q + 2Q) + fo 



(A.11) 



which yields the spectrum of (A. 4) as 



Ep±(q) = 2 - 5 - /„ - cos q - cos (q + 20) 
±^(cos q - cos {q + 2Q)Y + fl 



(A.12) 



The phase mode at £■ = corresponds to the eigenvector (1, 1) of the matrix 
(A.ll) at g = — Q, which imphes 



Ep_{-Q) = 2 - (5 - /o - /i - 2 cosg = 0, 



(A. 13) 



or 

5 = 5o{Q) -fo- fl. (A.14) 

Then, we have 

^P+(-g) = 2/i > 0, (A.15) 

Es+{-Q) = -2(/o - 2/i) = -2a, (A.16) 

and 

E5_(-g) = -2(/o + /i)<0. (A.17) 

For an analytic incommensurate SW, the sliding mode at = corresponds 
to the eigenvector (1, —1) of the matrix (A. 9) at g = —Q. Then, vanishing Eq. 
(A.16) yields a = so that fo = 2/i, and Eq. (A. 13) then yields 2cosg = 
2-S- 3/o/2. 



Otherwise, as explained in Sec. A. 2, for commensurate SWs represented by 
hyperbohc periodic cycles E — does not belong to the spectrum of the 
matrix (A. 9) for any q close to —Q, while for elliptic cycles = is in the 
spectrum of (A. 9) for some values of q close to (but away from) —Q {as E = 
is not at a band edge). In both cases, Eq. (A.16) yields that a ^ 0. Defining 
iasq = ^ — Q,^ small and real, it follows from Eq. (A. 10) that ii^ = is in 
the spectrum of (A. 9) if and only if the equation 
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Es+{^ - Q)Es-{^ - Q) = 4 cos' e - 4(2 - 5 - 3/o) cos Q cos ^ 

+ (2 - 5 - 3/o)' - 9f^ + 4cos2 Q - 4 = (A. 18) 

has a solution for ^ real. With cosQ taken from Eq. (A. 13), this equation 
becomes using (A. 6) 

(1 - cosO (2 + 2cos^ -{2-5- 3/o)(2 -5- 3/o/2 + a/2)) 

= a(3/o-«). (A.19) 

At leading order for ^ small and /o small, Eq. (A.19) yields 

^e'(4-(2-5)') = «(3/o-a). (A.20) 

Since < 5 < 4 for all Q 7^ if /o is small enough, there is a real solution 
for small ^ only when a > 0. Consequently, as a perturbative result valid 
at small amplitude, E = belongs to the spectrum of (A. 9) for some real 
1^1 > when a > 0, and does not when a < 0. Thus, the commensurate SW 
is represented by an elliptic periodic cycle when /o > 2/i and a hyperbolic 
periodic cycle when /o < 2/i. 



A. 3. 2 Band intersection with E = 
The determinant of the matrix (A. 7) is 

D{q- u;) = u''- {a\q) + a\q + 2Q) - 2/' + 6/1') 
+ ((a(g) - /o)(a(g + 2Q) - f,) - 9f^) x 

((a(g) + /o)(a(g + 2Q) + /o) - f^) . (A.21) 

Vanishing D{q; u) yields the intersections with the axis E = of the perturbed 
bands at wave vector q close to —Q. Thus, stability requires the zeros ujs{q) of 
D{q;uj) to be real for all wave vectors q close to —Q. At q = —Q Eq. (A. 13) 
yields a{—Q) = a{Q) = — (/o — /i), so that D{—Q;u) has a degenerate zero 
at a; = 0, while the second pair of zeros is given by 

cu',i-Q) = -4/i(/o - 2/1) = -4/ia. (A.22) 



Eq. (A.19) may also have solutions for ^ 2Q (g f» Q), but these should be 
discarded since we have assumed q ^ —Q in the perturbation calculation leading to 
Eq. (A.IO). 
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Thus, they are real only when a < 0, so that for small-amplitude SWs repre- 
sented by elliptic periodic cycles there is always a non-oscillatory instability 
(with purely imaginary a;^) at wave vector q = —Q. 

When q — —Q + ^ with ^ small, Eq. (A. 8) becomes using again (A. 13): 

a{TQ + ~ «(<3) T 2^ sin g = -(/^ + a) 2^ sin Q, (A.23) 

which yields for the determinant (A. 21) when also /i is small and a -C fi- 

D{-Q + 00)^00^- 4(2^2 sin^ Q - fia)uj'^ - 48C^a/i sin^ Q. (A.24) 

Thus, at order Eq. (A.21) has no real zeros when 

/la (/la + 8f sin^ Q) < 0. (A.25) 

For commensurate SWs represented by hyperbolic periodic cycles, a < 0. 
Then this condition becomes 

t2 ^ /l II 

^ 8sin2g' 

which proves that at small amplitude, the commensurate SWs represented by 
hyperbolic periodic cycles are unstable by oscillatory modulational instabilities 
with wave vectors q close to but not equal to iQ. 

Note that this proof does not hold for analytic incommensurate SWs because 
then q; = at order /q. In that case, the first-order instabilities found in Sec. 
4.2 for < Q < 7r/2 may occur outside the immediate neighborhood of a; = 
where the condition q ^ —Q assumed here is not valid. 
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